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INTRODUCTION 


Miniature  cryogenic  refrigerators  are  used  in  a 
variety  of  military  and  civilian  applications,  probably 
chiefly  to  cool  infrared  detector  elements  to  the 
appropriate  operating  temperature.  Long  wavelength 
sensors  operating  in  the  8  -  l^U/wvregion  commonly  re¬ 
quire  refrigeration  down  to  temperatures  in  the  4  -  30°K 
range ,  and  mid-wavelength  detectors  operating  in  the 
2-8 ^A/^region  usually  perform  satisfactorily  at  a 
temperature  in  the  70  -  8o°K  neighborhood. 

Historically,  the  70  -  80°K  requirement  was  met  by 
employing  open  cycle  liquid  nitrogen  systems  to  provide 
the  detector  refrigeration.  At  least  as  early  as  about 
I960  mechanical  refrigerators  were  successfully  used  in 
this  application.  Early  refrigerators  employed  were 
based  on  the  Joule -Thompson  process  in  which  high  pressure 
gas  was  expanded  through  a  nozzle  to  produce  either  a 
liquid  or  a  cold  vapor  capable  of  removing  heat  from  the 
detector.  Closed  cycle  Joule-Thompson  devices  were 
plagued  with  reliability  problems  stemming  both  from  the 
high  compression  loads  in  the  miniaturized  compressors 
and  from  the  stringent  contamination  requirements  for 
the  gas  in  order  to  avoid  blocking  the  nozzle. 

To  improve  the  reliability  situation,  classes  of 
refrigerators  which  employed  a  periodic  flow  heat  exchanger, 
called  a  regenerator,  and  a  piston-type  expansion  process 
were  developed.  Some  of  the  earliest  work  on  these 
regenerator-based  machines  was  reported  in  the  1940s,  and 
by  i960,  Stirling,  cycle  versions  were  already  at  a  reason¬ 
ably  mature  stage  of  development.  Over  the  next  decade 
nearly  a  half  dozen  variations  of  the  Stirling  cycle  were 
explored,  some  rejected  as  mere  laboratory  curiosities 
and  others  developed  to  a  high  level  of  sophistication.  All 
of  these  devices,  in  the  final  analysis,  were  limited  in 
their  capabilities  by  the  performance  of  the  periodic  flow 
regenerator. 


In  regenerates  heat  exchanger  matrix  fabrication  the 
search  for  suitable  materials  always  came  back  to  lead 
as  the  leading  candidate  since  it  (lead)  has  a  higher 
specific  heat  at  low  temperatures  than  most  other  materials. 
Even  with  lead,  however ,  the  rapid  decrease  in  specific 
heat  at  temperatures  below  about  15°K  severely  penalized 
refrigerator  performance  at  the  low  end  of  the  temperature 
range.  Careful  optimization  of  the  various  machines  re¬ 
sulted  in  systematically  reducing  the  "dead  end"  temperature 
over  a  period  of  years  to  values  below  10 °K. 

At  this  low  end  of  temperature,  however,  the  machines 
normally  provided  very  little  excess  capacity  to  account 
for  variation  in  performance  due  to  seal  degradation, 
slight  gas  contamination,  variation  in  crankcase  temper¬ 
ature,  or  any  of  the  multitude  of  off-design  conditions 
that  can  be  encountered  in  actual  user  environments. 
Consequently,  variations  in  machine  performance  would 
directly  translate  into  variations  in  the  temperature 
platform  provided  for  the  infrared  detector  element. 

Variations  in  the  detector  temperature  are  highly 
undesirable,  since  they  can  often  directly  show  up  as 
variations  in  the  output  signal.  Output  signal  variations 
can  cause  confusion  either  by  being  interpreted  as  valid 
information  to  be  decoded  or  by  masking  low  level  informa¬ 
tion  content. 

To  minimize  this  temperature  variation  problem,  the 
impetus  for  better  regenerator  design  and  the  search  for 
high  specific  heat  materials  continued.  Various  specific 
heat  anomalies  yielding  high  heat  capacity  values  (but  over 
very  narrow  temperature  bands)  have  been  identified.  Many 
of  these  have  been  tried  in  regenerator  experiments, 
usually  with  only  limited  success.  One  specific  heat 
effect  which  has  not  yet  been  tried  experimentally  is 
achieved  when  a  bulk  material  is  subdivided  into  an  ultra- 
fine  powder  form. 

This  specific  heat  enhancement,  which  is  associated 
with  the  inherent  high  surface  energy  of  powdered  material, 


has  been  reported  several  times  in  the  literature.  To 
this  point  in  time,  only  a  minimum  amount  of  investigation 
into  predicting  the  magnitude  of  the  effect  has  taken 
place  and  no  applications  have  been  attempted.  This 
surface  enhanced  specific  heat  effect  is  particularly 
attractive  for  refrigerator  application  since  it  is  not 
limited  to  a  narrow  temperature  band.  So  far  as  is  known, 
it  extends  from  the  lowest  temperatures  of  interest  for 
these  applications  up  to  probably  30°K  (or  warmer), 
depending  on  the  material. 

The  purpose  of  the  research  program  reported  herein 
has  been  to  determine  the  suitability  of  this  surface 
enhanced  specific  heat  for  use  in  low  temperature  regen¬ 
erator  applications.  Objectives  established  to  satisfy 
program  requirements  included: 

o  Develop  a  predictive  model  for  surface  enhanced 
specific  heat 

o  Search  literature  and  catalog  experimental  program 

results  of  measurements  of  this  surface  specific 
heat  effect 

o  Develop  a  model  for  predicting  performance  of 
a  regenerator  fabricated  using  finely  divided 
powders 

o  Apply  this  new  regenerator  model  to  predict 

performance  capability  of  a  refrigerator  using 
the  surface  enhancement  effect 

o  Catalog  current  state  of  the  art  in  small 
particle  production 

o  Develop  concepts  for  producing  practical 
regenerators  which  take  advantage  of  the 
surface  enhancement  effect 

o  Provide  recommendations  for  further  work  and 
a  system-oriented  test  plan  for  evaluating 
refrigerator  improvements  from  this  effect. 

Each  of  these  objectives  is  addressed  in  detail  in  the 
sections  following. 
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II 


SPECIFIC  HEAT  THEORY 


2 . 1  Background 

From  thermodynamic  considerations,  the  concept  of 
heat  capacity  is  best  described  as  that  amount  of  heat 
necessary  to  be  added  to  a  substance  in  order  to  raise 
its  temperature  by  one  measurement  unit.  Historical 
experiments  to  measure  the  specific  heat  of  materials 
at  moderate  temperatures  showed  the  surprising  result 
that  the  constant  volume  specific  heat  of  all  monatomic 
solids  (C„)  approaches  a  definite  values 

Cv  -  3  -  3ft.  =  (o  cj\w/kok«/*c  (i) 

where  ^5© is  Avogadros  number 

is  "the  Boltzman  constant 

and  Ti  is  the  gas  constant. 

This  relationship,  known  as  the  law  of  Dulong  and 
Petit,  has  been  found  to  hold  over  a  broad  temperature 
range.  Its  subsequent  theoretical  derivation,  which 
was  accomplished  utilizing  classical  statistical  mechanics 
principles,  was  considered  to  firmly  establish  it  as  a 
fundamental  law  of  nature.  The  relationship  was  further 
strengthened  when  experimental  results,  in  accord  with  the 

predictions  of  classical  theory,  demonstrated  that  diatomic 
solids  had  a  value  approximately  double  that  of  Eq.  1 
triatomic  solids  triple  ,  and  so  on. 

It  was  a  matter  of  considerable  concern,  then,  that 
subsequent  experiments  indicated  that  the  specific  heat 
of  solids  in  the  cryogenic  temperature  region  was  a  strong 
function  of  temperature,  in  sharp  disagreement  with  the 
predictions  of  classical  physics.  A  satisfactory  explan¬ 
ation  of  this  low  temperature  "anomaly"  awaited  the 
development  of  quantum  mechanics  concepts  early  in  the 
twentieth  century. 

»  defined  as  hear  capacity  per  unit  mass 
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The  turn  of  the  twentieth  century  marked  profound  advances 
in  our  understanding  of  science.  Planck's  brilliant 

empirical  approach  quantizing  the  energy  associated  with 
black  body  radiation  formed  the  basis  for  what  was  to 
become  the  quantum  theory  of  physics.  In  1906  Einstein 
applied  these  same  quantum  concepts  to  the  model  of 
specific  heat  with  enough  success  to  pave  the  way  for 
later  work  which  finally  explained  those  early  anomalous 
specific  heat  measurements.  The  most  influential  of 
the  specific  heat  models,  that  developed  by  Debye, 
accurately  predicted  a  low  temperature T* variation  in 
specific  heat. 

Before  developing  the  concept  of  improving  the  low 
temperature  specific  heat  of  a  material  by  finely  sub¬ 
dividing  it  into  an  ultrafine  powder  form  --  the  primary 
topic  of  this  study  program  --  it  is  useful  to  review 
the  Debye  specific  heat  analysis  for  bulk  material. 


Bulk  lattice  Specific  Heat 

The  first  law  of  thermodynamics  is  written  as 


Sq  =  «AU-v1?<JV 


where Sqis  heat  absorbed  by  a  system 

is  the  change  in  internal  energy 
"%Vis  the  work  done  by  the  system 
and  the  notation  SQ  indicates  the  quantity  is  an 
inexact  differential.  Expressed  as  a  function  of 
volume  (V)  and  temperature  CD. 
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Consequently , 

t 

Sq  -  (%\  iT-*-  IjCltfX  Jr'?l  <1  V .  ( 4 ) 

We  then  define  specific  heat  at  constant  volume  as: 


Eq.  5  will  be  the  basis  for  computing  specific  heat 
once  we  derive  an  expression  for  the  internal  energy  of 
a  solid  from  quantum  mechanics  principles. 

Quantum  theories  of  specific  heat  all  have  in  common 
the  postulate  that  the  total  internal  (vibrational) 
energy  of  a  crystal  of  NJ  lattice  sites  is  equivalent  to 
the  energy  of  a  system  of  3fs)  harmonic  oscillators. 
Consequently,  in  a  three  dimensional  solid  each  atom 
can  contribute  three  normal  modes  of  oscillation  to  the 
overall  vibrational  frequency  spectrum  which  we  shall 
show  determines  the  internal  energy  of  the  crystal. 

It  is  not  intuitively  obvious  that  so-called  vibration¬ 
al  modes  of  a  crystal,  which  most  readily  bring  to  mind 
a  wave  traveling  down  a  taut  string,  can  be  equated  to 
the  vibrations  of  a  three  dimensional  harmonic  oscillator, 
as  the  postulate  requires.  The  text  by  Dekker1  elegantly 
derives  the  mathematical  correspondence  between  the  two 
in  a  proof  which  won’t  be  repeated  here.. 

The  specific  heat  model  in  general  use  today  is  the 
Debye  equation  which  was  derived  only  a  few  years  after 
Einstein's  classic  work.  This  model  provides  highly 
accurate  predictions  of  specific  heat,  even  though  it 
was  derived  using  several  simplifying  assumptions,  result¬ 
ing  in  some  theoretical  shortcomings. 

Although  the  Debye  algorithm  can  be  derived  from 
ouantum  statistical  mechanics  principles,  better  insight 
to  the  physics  of  the  specific  heat  mechanism  is  gained 
by  the  following  less  rigorous  approach. 


/..>* 
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The  three  dimensional  time  dependent  wave  equation 
applicable  to  a  harmonic  oscillator  is, 

as  _  j_  sia 

"  c*  («) 

where  is  the  instantaneous  wave 

deflection  at  a  location 
O^Hj-^in  space,  "t  is  time, 

and  C*  is  the  propagation  velocity 

in  the  medium. 

A  solution  to  this  equation  is  the  wave  of  frequency 
•V  described  by: 

/^s,itk^n>4Trv/L^si>(mY‘n‘  y/1)* 
sm  ^  Cos  2n  V^:  ^  ^ 7  * 

provided  that  /*'-».>  /"*«¥  and  are  all  positive 
integers  ^ l ,  and  the  solid  is  in  the  form  of  a  cube 
of  side  L . 

The  frequencies  allowed  by  Eq.  7  take  the 
form: 


This  relationship  exhibits  the  fact  that  only  integral 
numbers  of  half-wave  lengths  are  allowed  to  propagate 
through  the  crystal  lattice. 

In  order  to  compute  the  internal  energy  from  which 
to  derive  the  specific  heat,  we  are  interested  in 
determining  the  possible  modes  of  vibration  in  the 
cube,  multiplying  each  by  the  harmonic  oscillator 
energy  associated  with  it,  and  then  summing  all  such 
mode  energies  to  determine  the  total  for  the  crystal. 

Since  the /A^,  of  Eq.  6  are  required  to  be  all 
positive,  we  can  associate  them  in  integer  space  with 


-  - 


i 

•  Jtii 
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the  first  octant  of  a  sphere  of  radius^,  where 


*  /*£  -4/v^  .  (9) 

Now  since  each  positive  integer  value  gives  rise  to  a 
discrete  mode  of  vibration,  the  number  of  such  modes 
in  frequency  space  lying  between  H  and  designated 

dtaw*  can  be  approximated  by  the  volume  of  such  an 
octant.  Consequently, 

*£*§[■)  vV?  do, 

where  V ,  the  volume  of  the  cube,  is  L?  . 

The  approximation  sign  in  this  equation  is  only  to  alert 
the  reader  that  the  number  of  integer  points  in  such  an 
octant  is  not  rigorously  determined  by  such  an  integration 
technique.  This  concept  will  be  pursued  in  some  depth 
later. 

It  must  also  be  pointed  out  that  elastic  waves  in  a 
crystal  propagate  either  as  longitudinal  or  as  transverse 
waves,  and  different  velocities  Cr  3  CL  may  be  associated 
with  the  two.  Since  there  are  two  orthogonal  transverse 
waves  allowed,  we  modify  Eq.  10  to  the  more  general  form 
(without  carrying  forward  any  longer  the  approximation 
symbol) : 

*  4*  V  ^7  ^  -V*  Jv.  ( 11) 

We  see  immediately,  from  our  earlier  assumptions, 
that 

% 

\  -  3  Isl  (12) 

whereVpis  the  so-called  Debye  frequency,  which  cuts  off 
the  number  of  modes  added  up  by  this  integral  at  a 
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frequency  so  the  sum  equates  to  the  3N  allowable 

modes  in  the  crystal.  _ _ 

Since  Planck  showed  that  the  average  energy of 
a  harmonic  oscillator  at  temperature  l  is: 


_  _ _ 

% « -) 


(13) 


wherein,  is  the  Planck  constant, 
then  the  total  internal  energy  of  the  crystal  is  deter¬ 
mined  conceptually  by  multiplying  each  "shell"  of  modes 
described  by  Eq.  11  by  the  Planck  relation  applying 
to  that  shell,  and  adding  them  all  up.  Thus , internal 
energy  of  the  crystal  is  expressed  mathematically  as. 


u  =  c  ^ .  uv 

Wo 

If  we  define  the  Debye  temperature,  e„  ,  as  ^ 

then  we  find  from  Eq.  1^  that  at  low  temperatures 
the  internal  energy  per  mole  is: 


,  ,  ,  i  *» 
U  -  5-TT 


6>r 


(15) 


and  from  Eq.  (5) 


(16) 


Eq.  16  is  the  desired  relationship  for  bulk 
lattice  specific  heat  at  low  temperature,  and  shows  the 
T  -cubed  dependence  noted  earlier. 

In  the  Debye  approach  just  outlined,  we  have  implicitly 
assumed  that  the  crystal  could  be  represented  as  a  con¬ 
tinuum,  and  that  all  waves  propagate  as  though  they  were 
embedded  in  a  perfectly  elastic  medium.  Although  this 
car:  be  shown  to  be  a  good  approximation  for  low  frequencies 
where  the  wavelength  of  such  a  disturbance  is  long  compared 


with  the  interatomic  spacing  of  the  crystal,  for 
shorter  wavelengths  comparable  to  the  lattice  constant 
of  the  material  such  a  model  breaks  down. 

Consequently,  the  approach  employed  in  arriving  at 
Eq.  14  is  generally  valid,  but  the  representation 
of  bU-y;  as  a  continuous,  simply  derived  function  is  mis¬ 
leading.  If  the  lattice ,  instead  of  being  treated  as  an 
elastic  continuum  is  considered  to  be  made  up  of  an 
infinite  array  of  mass  points  which  interact  with  one 
another  via  harmonic  restoring  force,  and  it  is  further 
recognized  that  although  nearest  neighbor  interactions 
predominate  some  coupling  to  further-removed  lattice 
sites  may  exist,  then  a  relationship  for  )  cannot  in 
general  be  derived  in  closed  form. 

Even  where  appropriate  simplifying  assumptions  are  applied 

to  the  general  discrete-lattice  problem  it  is  found  that 

the  velocity  of  propagation  varies  with  the  wavelength  of 

the  disturbance  (in  contrast  to  the  constant  value  for  Cy 

andC^used  in  the  Debye  theory).  Further  it  is  found  that  the 

allowable  range  of  frequencies  which  can  propagate  through 

the  crystal  is  separated  into  a  low  frequency  band 

(coinciding  with  the  so-called  acoustical  waves  of  the 

Debye  theory)  and  a  high  frequency  band  (associated  with 

wavelengths  comparable  to  the  infrared  energy  spectrum) . 

These  are  referred  to  respectively  as  the  "acoustical" 

and  "optical"  branches  of  frequencies.  Since  an  energy 

gap  (or,  more  accurately,  a  band  of  "forbidden"  frequencies) 

exists  between  the  acoustical  and  the  optical  branches, 

it  is  apparent  that  a  plot  of  against  V  for  real 

substances  will  show  a  discontinuity  in  the  gap  region. 

Figure  1  contrasts  the  simple  continuum  assumption  of 

Debye  with  the  approximate  model  for  a  simple  ionic 

2 

crystal  like  sodium  chloride  . 

Although  rigorous  solutions  to  the  form  of 
cannot  be  obtained,  in  general,  the  lattice  dynamics 
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community  has  made  considerable  progress  in  the 
analytical  approach  to  modeling  its  behavior  for  some 
of  the  simpler  crystal  structures. 

2 . 3  Surface  Contributions  to  Specific  Heat 

The  Debye  formulation,  although  an  oversimplification 
of  the  physics  of  the  situation,  represents  an  outstanding 
success  in  applying  the  then-new  concepts  of  quantum 
mechanics  to  a  problem  which  couldn't  be  solved  with 
classical  techniques.  That  model,  as  just  outlined, 
was  derived  from  considerations  of  the  three-dimensional 
characteristics  of  an  elastic  continuum  as  an  approximation 
to  a  crystal  lattice. 

Since  real  substances  exhibit  a  two-dimensional  surface 
and  not  solely  an  infinite  three-dimensional  lattice,  we 
may  well  inquire  as  to  whether  or  not  there  can  be  a 
comparable  energy  contribution  to  the  overall  internal 
energy  stemming  from  the  crystal  surface  in  addition  to 
the  contribution  from  the  bulk  lattice.  Anticipating  the 
discussion  to  follow,  we  will  find  that  in  exact  analogy 
to  Debye's  derivation  a  surface  contribution  to  internal 
energy  (and  thus  via  Eq.  5  to  total  specific  heat) 
can  be  predicted.  As  we  would  intuitively  expect  from 
an  analogy  to  Eq.  10  }  the  contribution  turns  out  to  be 
proportional  to  the  surface  area  of  the  sample  and  is 
only  significant  under  conditions  of  very  large  surface- 
to-volume  ratios  and  at  low  temperature. 

It  will  be  remembered  that  Eq.  10  was  obtained  by 
estimating  the  number  of  volume  units  in  the  first 
quadrant  of  a  sphere  in 'V- space.  The  estimation  method 
used  is  accurate  enough  for  large  volumes,  but,  as  noted, 
is  not  exactly  equal  to  the  actual  count  of  non-negative 
integer  points  in  the  quadrant  when  the  surface-to-volume 
ratio  of  the  substance  is  large.  Consequently,  the  actual 
number  of  allowable  modes  is  greater  than  that  indicated 
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by  this  equation.  This  can  be  conveniently  illustrated 
with  a  two-dimensional  example: 


AREA  APPROXIMATION  TO  INTEGER  POINTS 


♦  •  •  • 


•  •  • 


•  • 


•  • 


4 


Number  of  integer  points  =  20 
Area  of  figure  =  12  units 

Half  perimeter  of  figure  =  7  units. 


FIGURE  2 

If  the  area  of  Figure  2  is  used  to  approximate 
the  number  of  non-negative  integer  coordinate  points  of 
the  figure,  the  result  is  in  error.  Approximately  one- 
half  the  points  on  the  perimeter  must  be  added  to  this 
two-dimensional  area  value  as  a  correction. 

In  the  three-dimensional  case,  the  volume  estimation 
is  similarly  in  error,  and  half  the  surface  area  must  be 
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added  as  a  correction.  E.  W.  MontrollJ  was  one  of  the 

first  workers  to  recognize  this  discrepancy  in  Debye’s 

specific  heat  model.  Drawing  on  earlier  work  from  the 

4 

field  of  acoustics  ,  he  derived  a  frequency  distribution 
function  which  accounted  for  lattice  contributions  plus 
the  correction  for  allowable  vibrational  modes  associated 
with  the  large  surface  area  of  powders  and  related 
geometries.  Montroll's  resulting  form  of  Eq.  10  is, 


. . L  +  negligible  edge 

4ttv  Q  **•  XTT  oCj,  correction  terms. 


(16) 


where  S  is  the  total  surface  area  of  the  substance. 
Using  the  same  Debye  approach  illustrated  in  Eq.  14  we 
find 


K)  -to. 


(1?) 


Jl3) 


^  -A 

where  a  Cu  ■*  “2.  CT 

is  the  Riemann  zeta 
value  of  3- 


(18) 


function  evaluated  for  a 


For  a  mole  of  the  substance  this  can  be  expressed  as; 


(19) 


T  3 


where  the  l  coefficient  has  the  bulk  Debye 
value  of; 

„  /  / 

p  9  gg  CAu/houE/jc  • 

.  .  -3  _z 

Ratiomg  the  (  term  to  the  l  term  we  finds 

g&fti*  _ _ 


(20) 


(21) 
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or 


(21a) 


where 


Thus,  the  first  Montroll  coefficient  is  found  to  be: 


.on  p*  cl ) 


(22) 


where  M  is  the  atomic  weight 

and  ^  is  the  bulk  density  of  the  material. 

If  we  accept  Montroll's  statement  that  it  can  be  shown 
that  <  CSL  4  i 


then 


If  specific  surface  area(5*)  is  defined  as: 

.  s.  -  Su  =,  ±.  ) 

S  ~  AA.  ~  |»V  f  L  V  / 


(22a) 


(23) 


where  fit V#  is  the  mass  of  the  sample  , 
then  an  alternate  equation  for  this  Montroll  coefficient 
is : 


(24) 


For  discussion  purposes  this  can  be  rearranged  as: 


c<  «C 


(25) 


It  is  not  surprising  to  find  that  Eq.  25  indicates 
that  potentially  the  largest  coefficient  value  for  the 
surface  contribution  to  specific  heat  will  be  found  for 
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heavy,  dense  materials  having  the  lowest  possible  Debye 
temperature.  The  Debye  temperature  is  obviously  the 
most  important  selection  parameter,  according  to  Eq.  25  » 
and  we  see  that  bulk  density  is  relatively  more  important 
than  atomic  number  in  choosing  a  material  expected  to  show 
large  surface  specific  heat  effects  in  the  finely  sub¬ 
divided  state. 

It  is  appropriate  to  briefly  summarize  the  preceding 
mathematical  derivations.  If  we  consider  a  continuous 
three-dimensional  lattice  structure,  a  standard  Debye- 
type  derivation  from  quantum  mechanics  principles  shows 
that  the  specific  heat  of  solids  varies  with  the  cube  of 
temperature  in  the  cryogenic  region.  If,  however,  we 
consider  a  special  configuration  of  the  solid  wherein  a 
large  surface-to-volume  ratio  exists  --  e.g.,for  finely 
subdivided  powders  or  thin  films  --  the  Montroll  approach 
predicts  an  additional  contribution  to  specific  heat  which 
varies  as  the  square  of  temperature;  consequently,  this 
surface  contribution  to  specific  heat  does  not  decrease 
at  low  temperatures  as  rapidly  as  the  Debye  lattice  term. 

The  magnitude  of  this  surface  contribution  can  be  pre¬ 
dicted  from  physical  characteristics  of  the  substance, 
as  indicated  by  Eq.  25  . 

Montroll's  approach  is  subject  to  the  same  theoretical 
criticism  as  was  levied  on  the  simple  Debye  model  discussed 
earlier;  e.g.,  an  elastic  continuum  must  be  assumed  instead 
of  a  discrete  lattice  of  coupled  mass  points,  and  simple 
periodic  boundary  conditions  are  assumed  in  order  to  allow 
solution  of  the  wave  equation.  In  addition,  it  has  been 
noted  that  this  Montroll  approach  does  not  take  into  account 
the  fact  that  longitudinal  mode  waves  and  transverse  mode 
waves  may  undergo  a  mode  interchange  upon  reflection  at 
crystal  boundaries,  and  that  under  certain  conditions 
additional  surface-localized  and  edge-localized  modes  are 
allowable . 


The  net  result  is  chat  the  Montroll  approach  also 
under-counts  the  number  of  allowable  vibration  modes 
(although  not  as  flagrantly  as  the  Debye  model),  and 
should  tend  to  err  in  the  direction  of  predicting  too 
low  a  value  for  specific  heat.  This  has  been  verified 
in  the  few  reported  measurements  on  finely  divided  powders 
which  can  be  correlated  to  specific  heat  predictions  of 
the  model^. 

There  have  been  many  other  approaches  to  modeling 
surface  specific  heat  in  the  twenty-odd  years  since  the 
work  of  Montroll.  At  present  none  of  the  lattice  dynamics 
models  developed  in  the  last  decade  appear  to  be  superior 
to  the  simple  model  of  Eq.  25  in  predictive  ability. 

It  is  recognized,  however,  that  lattice  dynamics  approaches 
will  eventually  develop  an  accurate  algorithm  for  comput¬ 
ing  surface  specific  heat  contributions.  The  solution  to 
the  problem  may,  in  fact,  already  be  at  hand,  but  unveri¬ 
fied  due  to  a  lack  of  reported  measurements  on  fine  powders; 
however,  any  equations  stemming  from  lattice  dynamics  work 
are  likely  to  be  complex  and  tedious  to  use. 


Because  of  the  fact  that  lattice  dynamics  approaches 
do  not  presently  have  a  clear  advantage  in  predictive 
ability  over  the  simple  Montroll  model,  we  will  not  dwell 
on  that  approach  in  any  depth.  It  would  be  remiss,  however, 
to  gloss  over  the  subject  without  illustrating  a  recent 
effort  to  develop  a  simplified  lattice  dynamics  equation 
for  the  surface  specific  heat.  Naugle  and  Allen^  have 
applied  their  approximate  algorithm, 


=  Tffc  L  ~  1A  T.-n-fc.T  /  _ 


AjAeawhc  ^ 


(26) 
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*  -mV  . 


where 


-fby  =  w  (*4) 


, .  1*. 

0 


and  LO^  =  'Z'K'V^  5 

to  noble  gas  solids  and  alkali  halides  with  the  result 
they  have  closely  matched  values  computed  from  rigorous 
lattice  dynamics  equations.  The  Wo.  L  i  anc^ 

cO^T^in  this  equation  correspond,  respectively,  to  surface 
vibration  frequencies  of  the  acoustical  and  the  optical 
branches  after  the  formation  of  a  surface,  and  before  the 
surface  is  created  (superscript  b)  --  i.e.,  in  the  bulk. 

From  a  theoretical  standpoint  this  mathematical 
approach  is  very  satisfying.  From  an  experimental  stand¬ 
point,  however,  this  equation  and  its  more  rigorous 
kindred  are  no  better  at  predicting  surface  specific  heat 
than  the  simple  model  outlined  earlier. 

The  equation  resulting  from  the  Montroll  approach  is 
convenient  to  use,  but  as  might  be  expected,  it  has  some 
shortcomings.  First,  it  is  not  really  the  most  accurate 
predictor  available.  Based  on  the  few  published  data  sets 
to  date ,  the  most  accurate  algorithm  is  one  derived  by 
two  different  workers  independently,  circa  I960'7.  In  this 
more  accurate  Depuis-Mazo-Onsager  (DMO)  model,  the  same 
Debye  approach  --  involving  assumptions  of  an  isotropic 
elastic  continuum  --  was  employed.  The  primary  difference 
from  the  Montroll  work  was  that  the  DMO  model  attempts 
to  impose  more  realistic  boundary  conditions  on  the  wave 
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equation,  and  it  recognizes  that  introduction  of  a 
surface  into  the  continuum  "scrambles"  wave  reflections 
in  such  a  way  that  resulting  vibrations  can't  be  analyzed 
in  terms  of  pure  transverse  and  longitudinal  modes. 

The  form  of  this  equation  for  surface  specific  heat  Gy 

s  ./ ,  STi 

r  qf*  y3j  -■*  -i  /  .*»  -  ■»  '  ‘  (??) 


C?  cl  (cl  -  CTl  ) 


where  the  symbols  are  as  earlier  defined. 

In  order  to  employ  the  DMO  algorithm,  it  is  necessary 
to  know  the  low  temperature  elastic  constants  of  the 
material  of  interest.  Although  CT  (average)  can 
apparently  be  estimated  from  the  experimental  Debye 
temperature ,  ©&  ,  in  the  temperature  range  of  interest 
Ck.  is  not  generally  available.  In  addition,  the  low 
temperature  ©p  varies  widely  as  I- — >  O  ,  so  even  the 
determination  of  Cx  from  this  relationship  is  less  than 
satisfactory. 

The  most  recent  analytical  treatment  of  surface  enhanced 

8 

specific  heat  seems  to  be  that  of  Baltes  and  Hilf  .  They 
developed  the  algorithm,  a 

^  Hit  n?  F*  (28) 

where  -r~ — -  »  ^  denotes  the  particle  radius,  \4w\ 

w  4tt  7L~ T  •* 

is  the  mole  volume  of  the  sample ,  and 
of  t  As  and  Si  are  associated  with 
spherical  Bessel  functions. 

Although  recent  measurements  show  that  this  algorithm 
is  potentially  very  accurate,  like  the  DMO  model  it  is 
very  difficult  to  use  because  the  effective  sound  velocity, 
Gx  ,  varies  strongly  at  low  temperature  and  can  best  be 
determined  by  measuring  the  variation  of  the  Debye  temp¬ 
erature  itself  in  the  range  of  interest.  Consequently, 
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the  equation  is  probably  not  very  useful  for  our  purposes 
as  a  predictor,  since  knowledge  of  the  Debye  temperature 
variation  is  not  usually  available  for  a  given  material. 
Consequently,  we  again  are  led  back  to  the  simpler 
Montroll  model. 

o 

In  a  recent  survey  of  the  ability  of  the  various 
models  to  accurately  estimate  the  surface  specific  heat 
value  for  dielectric  powders,  it  was  noted  that  results 
from  the  Montroll  approach  underestimate  the  surface  con¬ 
tribution  to  specific  heat  (i.e.  give  a  value  too  low)  by 
a  factor  of  4  -  6.  The  more  complicated  DMO  model  can 
improve  upon  the  Montroll  accuracy  by  only  about  45  percent, 
and  at  this  point  in  time  there  are  no  other  models  (in¬ 
cluding  the  most  sophisticated  lattice  dynamics  algorithms) 
which  can  do  any  better  without  requiring  extensive  know¬ 
ledge  of  low  temperature  crystal  parameters. 

2 .4  Electron  System  Contribution  To  Surface  Specific  Heat 

The  preceding  derivations  used  as  a  starting  point 
the  bulk  lattice  and  "bulk  surface"  contributions  to 
specific  heat.  It  must  be  recognized  that  other  mechanisms 
can  also  contribute  to  specific  heat  of  a  substance;  e.g. 
order-disorder  transitions  in  the  crystal  lattice,  order- 
disorder  transitions  in  the  magnetic  spin  system  of  the 
substance,  and  energy  level  changes  in  the  conduction 
band  electron  system,  to  name  the  most  important.  The 
latter  effect  is  only  significant  in  metals,  due  to  the 
broad  range  of  energy  states  available  in  the  conduction 
band . 

Since  electrons,  being  particles  of  half-integral  spin 
values,  are  represented  by  anti symmetrical  wave  functions, 
we  must  apply  Fermi-Dirac  statistics  in  order  to  compute 
average  values  for  the  conduction  band  energy  (and  con¬ 
sequently  the  electron  contribution  to  specific  heat) . 

Rather  than  repeat  the  derivation  for  specific  heat  of  an 
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electron  "gas"  as  found  in  the  literature,  the  main 

difference  from  the  bulk  lattice  case  can  be  illustrated 
as  follows. 

It  can  be  shown  from  quantum  mechanical  principles10 
that  m  analogy  to  Eq .  10  the  number  of  allowable  energy 
states  U((£)  for  electrons  in  the  range  between  £  and 

fc^dEis 


E  * 


4 7rV ^ 


(29) 


where  /wv^is  the  electron  mass 
and  it  is  recognized  that  each  electron  of 
opposite  spin  represents  a  distinct  energy 
state . 

The  Fermi-Dirac  distribution  law  states  that  the 
probability  of  finding  an  electron  in  any  one  of  these 
energy  states,  £  ,  is 


-f60  =  ,  (30) 

where  Ep  is  the  Fermi  energy  level  of  the 
crystal „ 

Consequently,  the  average  energy  per  mole  of  the  system 
can  be  shown  to  be: 


E.  '  ^  dc  (31) 

°r  ,  ,«•  Uir  E*  V 

’V\  SL(B'e”^[  dE  .  Ola) 


The  solution  to  this  integral  is 

E  *  E0  ■+  -^TTZ  ( 31b) 
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which  yields  the  specific  heat  equation  of  interest: 

/  3  k/v  L\^ 


^  3  3  A? 

commonly  written 


T 


(32) 

(33) 


Thus  we  conclude  that  the  specific  heat  of  the 
electron  system  at  low  temperature  is  proportional  to 
the  first  power  of  the  absolute  temperature.  It  is  known 
that  because  of  the  assumptions  inherent  in  the  model  the 
numerical  value  of  measured  results  sometimes  deviates 
markedly  from  the  computed  value.  However,  the  expected 
linear  proportionality  to  l  does  indeed  occur. 

Now,  if  we  apply  the  same  rationale  used  earlier,  the 
density  of  state  function  of  Eq.  29  has  again  under¬ 
counted  the  number  of  allowable  electron  quantum  states 
associated  with  a  finely  divided  powder  or  similar  sub¬ 
stance  showing  a  large  surf ace -to-volume  ratio.  Montroll 
showed  that  for  this  case  Eq.  29  must  be  rewritten  as 

9?d=)<i  e  -  [|r e''4-*  i  E  .  <34) 


Upon  performing  similar  steps  as  just  indicated  to 
compute  first  the  internal  energy  and  then  the  specific 
heat,  Montroll  showed: 
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It  will  be  noted  that  upon  rearrangement  the  first  term 
is  identical  to  the  expression  of  Eq.  32  .  In 
analogy  with  that  formula,  we  express  the  enhanced  specific 
heat  as: 
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and  where  is  the  average  number  of  conduction  electrons 
per  lattice  site. 

Numerical  calculations  show  that  this  6  term  is  small 
for  most  metals;  consequently  the  enhanced  electron  con¬ 
tribution  is  a  less  important  small  particle  effect  than 
the  surface  vibrations  described  by  the  oCT7"  term  predicted 
by  Montroll. 

2 . 5  Survey  of  Reported  Measurements 

The  interest  in  surface  contributions  to  specific  heat 
at  low  temperature  was  perhaps  stimulated  by  early  measure¬ 
ments  on  microcrystals  by  Giauque  and  Archibald11, 

reported  in  1937.  In  the  early  1950 ' s  a  flurry  of  theoret¬ 
ical  work  in  surface  thermodynamics,  highlighted  (at  least 
for  our  purposes)  by  the  models  of  Montroll, spurred  some 
additional  experimental  measurement  work  on  finely  divided 
powders . 

12 

Patterson,  Morrison,  and  Thompson  reported  in  1955 
that  they  were  unable  to  detect  the  *  Y2,  term  in  low 
temperature  measurements  on  ,  but  did  detect  this 

T  -squared  effect  on  measurements  of  finely  divided  koCl 
It  is  generally  conceded  now  that  the  effect  in  X icv 
was  undoubtedly  present,  but  masked  by  experimental  errors. 

In  the  case  of  tVJaCJL  ,  these  researchers  found  perhaps 
the  first  indication  that  the  magnitude  of  the*T*- squared 
term  measured  experimentally  in  dielectric  powders  is 
significantly  larger  than  the  value  predicted  by  Montroll's 
relationship  of  Eq.  24  .  Their  actual  reported  measure¬ 
ment  values  are  suspect,  however,  because  of  several 
experimental  difficulties  they  were  not  able  to  satisfactorily 
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resolve.  Some  significant  errors  were  made,  for  example, 
in  the  measurement  of  powder  surface  area.  The  difficulty- 
in  determining  the  area  of  powder  samples  will  be  discussed 
in  a  later  section  of  this  report,  but  it  is  appropriate 
to  introduce  the  methodology  at  this  point  since  it  bears 
on  the  problems  encountered  by  these  authors. 

In  1938,  Branauer ,  Emmett  and  Teller^  developed  a 
comprehensive  gas  adsorption  theory  which  provided  a  means 
of  determining  the  surface  area  of  finely  divided  solids 
from  isothermal  gas  adsorption  measurements.  The  technique 
involves  a  determination  of  the  mass  quantity  necessary  to 
coat  the  solid  with  a  monolayer  of  the  gas.  Even  with  the 
advent  of  X-ray  diffraction  methods  and  scanning  electron 
microscope  technology,  this  BET  adsorption  isotherm  method 
still  today  remains  the  standard  for  determining  the  surface 
area  of  fine  particles.  In  the  measurements  of  Patterson, 
et  al ,  however,  an  anomalous  behavior  of  the  adsorbent- 
adsorbate  system  resulted  in  a  non-linear  BET  plot  even 
in  the  range  of  gas  adsorbate  pressures  previously  believed 
suitable  for  measurements  of  this  type,  and  gave  rise  to 
the  error.  The  difficulty  was  subsequently  resolved  by 

,  .  ih. 

Maclver  and  Emmett  (of  the  original  BET  team)  and  apparently 
has  not  hampered  later  workers. 

Following  Patterson's  work,  Lien  and  Phillips^  reported 
the  first  data  on  small  particle  heat  capacity  which  can  be 
considered  generally  reliable.  Their  data  on  in  the 

1.5°K  to  4°K  range  clearly  showed  the  0<T  effect  predicted 
by  Montroll.  The  M^O  sample  was  in  the  form  of  ultra 
fine  particles,  with  a  cube  edge  of  100  K.  The  co¬ 
efficient  determined  by  them  was  larger  than  the  Montroll 
prediction  by  nearly  an  order  of  magnitude. 

-1  /r 

In  1966  Barkman,  Anderson,  and  Brackett±  performed 
a  carefully  designed  set  of  experiments  to  show  conclusively 
both  that  the  T-  squared  term  actually  exists  in  the  finely 
divided  samples  and  that  its  magnitude  is  proportional  to 
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surface  area.  Their  classic  work  with  k*.OL  ,  over  a 
broad  range  of  sample  surface  areas  and  for  a  temperature 
range  from  4.5°K  to  18.1°K,  demonstrated  the  existence  of 
the  surface  enhanced  specific  heat  and  further  showed  that 
the  DM0  equation  is  its  most  accurate  predictor,  with  the 
Montroll  equation  giving  values  of  the  surface  contribution 
too  low  by  a  factor  of  4  -  6, 

In  the  past  few  years  there  has  been  a  flurry  of 

interest  in  measurements  of  the  surface  specific  heat  of 

metal  particles.  In  a  very  excellent  experimental  program 

"  17 

m  Germany,  Comsa,  Heitkamp  and  Rade  measured  the  specific 

heat  of  palladium  ultrafines  of  30  and  66  &  diameters.  > 

Their  work  has  dramatic  implications  on  several  aspects 

of  this  regenerator  improvement  study,  and  will  be  cited 

in  several  following  sections. 

Their  work  on  metal  powders,  similar  to  the  Barkman 
work  on  dielectrics,  shows  conclusively  that  a  strong 
surface  enhancement  in  specific  heat  is  associated  with 
particle  size.  Comparable  to  the  earlier  indications, 

the  Montroll  equation  is  again  found  to  predict  values  j 

too  low  by  a  factor  in  the  4-6  range.  Their  results 

did  not  detect  that  the  electron  specific  heat  for  the  1 

ultrafines  was  enhanced  over  the  value  associated  with 
bulk  palladium.  This  is  not  inconsistent  with  the  analyti¬ 
cal  result  noted  earlier:  that  such  an  enhancement,  if  it 
exists,  would  be  numerically  small  in  comparison  with  the 
terms  attributed  to  surface  mode  vibration. 

Several  years  before  the  palladium  measurements  were 
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reported,  Novotny  and  Meincke  also  verified  the  surface 
enhancement  term  in  37  A  and  60&  particles  of  lead,  and 
22  &  particles  of  both  lead  and  indium.  Especially  in  the 
case  of  lead,  the  surface  enhancement  contribution  was 
considerably  smaller  than  we  would  expect  from  a  Montroll- 
type  analysis.  Consequently,  it  was  of  priority  interest 
to  determine  whether  the  measurements  indicate  a  valid 
deviation  from  the  algorithm  of  Eq.  22,  thereby 


25 


raising  doubts  as  to  the  generality  of  the  Montroli  con¬ 
cept,  or  whether  peculiarities  of  the  experimental  technique 
used  will  allow  us  to  accept  these  results  without  in¬ 
validating  the  Montroll  theory.  These  questions  are 
addressed  in  Section  2.7. 

As  can  be  inferred  from  the  brief  length  of  the  fore¬ 
going  summary,  there  have  been  very  few  reported  measurement 
programs  on  the  specific  heat  of  small  particles.  This  is 
not  due  to  a  lack  of  scientific  interest  in  surface  thermo¬ 
dynamics  of  powders  at  low  temperatures,  since  the 
literature  contains  many  recent  theoretical  treatments  of 
the  problem.  We  can  conclude  that  the  lack  of  data  is  due 
at  least  in  part  to  the  difficulty  in  performing  the  required 
measurements.  These  difficulties  can  be  grouped  into 
several  categories,  each  of  which  will  be  treated  in  more 
detail  later  in  this  report: 

1)  The  production  of  sub-micron  powders  is  difficult  at 
best,  and  is  not  currently  done  on  a  commercial  basis; 
consequently,  a  researcher  traditionally  had  to 
produce  his  own  powders  for  use  in  a  measurement 
program  of  this  nature. 

2)  All  powders  in  this  size  range,  and  particularly 
metals,  are  typically  very  reactive,  both  in  the 
sense  that  they  tend  to  agglomerate  into  larger 
particles  (+v'us  losing  the  desired  surface  character¬ 
istics)  and  that  they  can  be  pyrophoric,  tending 
toward  spontaneous  combustion  and  even  explosion. 

3)  The  inherent  difficulty  of  the  measurements  to  be 
made  including  in  particular  surface  area 
determination,  sample  temperature  determination, 
and  determination  of  effective  heat 

input  during  the  calorimetry  process  can  j eopardize 
even  a  well  designed  specific  heat  experiment  program. 
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2.6 


The  Question  of  Cluster  Structures 

There  are  many  different  physical  processes  which 

tend  to  produce  a  numerical  deviation  from  the  Debye  TJ 

specific  heat  law  at  low  temperatures.  Most  of  these, 

such  as  the  low  temperature  order-disorder  transitions 

and  the  various  low  temperature  magnetic  alignment  effects 

occur  only  in  a  relatively  small  class  of  materials  and 

only  over  a  very  narrow  temperature  range,  on  the  order 
o  19 

of  a  few  K  .  In  contrast,  the  surface  enhanced  excess 

specific  heat  which  is  the  primary  topic  of  this  study 

is  generally  applicable  to  all  substances  which  can 

be  subdivided  fine  enough  that  the  coefficient  of 

Eq.  22  is  measurable.  Additionally,  the  surface  effect 

2 

is  a  "broad  band"  enhancement,  varying  as  T  ,  and  is 

not  limited  to  a  narrow  temperature  range. 

Although  we  can  assume  that  a  broad  band  effect  is 

ultimately  more  useful  for  application  as  a  heat  exchanger 

(regenerator)  matrix  material  in  a  cryogenic  refrigerator, 

it  is  reasonable  to  suggest  that  selected  narrow  band 

materials  might  be  useful  in  the  cold  end  segment  of  the 

"segmented  regenerator  model"  introduced  later  in  this 

report.  Consequently,  some  effort  was  invested  during 

this  study  to  screen  other  specific  heat  anomalies  for 

their  potential  application  to  the  cold  segment.  In  the-.  ' 

process  of  that  survey,  an  interesting  coincidence  was 

discovered  which  may  allow  us  to  take  advantage  of  the 

Montroll-type  enhancement  without  the  difficulty  of 

having  to  work  with  finely  powdered  materials. 

2  0 

In  1961  Schroder  published  an  analytical  model  to 
describe  a  deviation  from  the  expected  linear  and  cubic 
specific  heat  terms  for  certain  metal  alloys  containing 
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small  ferro  magnetic  clusters  which  caused  the  entire 
lattice  to  exhibit  evidence  of  superparamagnet ism. 

Under  the  assumption  that  these  clusters  oscillate 
around  a  particular  preferred  orientation,  the  author 
showed  that  for  large  clusters  a  specific  heat  term  CtL 
of  the  Einstein  form, 

(38) 

should  be  found  experimentally.  In  this  equation, 
f\  equals  the  number  of  degrees  of  freedom 
for  the  cluster, 

Cg  is  the  Einstein  specific  heat  function 
and 

"Tg  is  the  Einstein  temperature  (defined  in  a 
manner  analogous  to  the  Debye  temperature). 

He  further  argued  that  for  these  ferromagnetic 
clusters  the  effective  Einstein  temperature  would  be  well 
below  1°K,  so  that  the  Einstein  function  would  reach  a 
constant  value  at  temperatures  above  2  or  3°K.  By  this 
analysis  Schroder  explained  an  anomaly  first  discovered 
in  a  series  of  Ni-Cu  alloys  wherein  experimental  specific 
heat  measurements  at  low  temperatures  were  apparently  fit 
by  an  equation  of  the  type , 

T5  .  .  (39) 

It  seemed  an  interesting  coincidence  that  some  ten  years 
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earlier  Jura  and  Pitzer  had  predicted  on  solely  theo¬ 
retical  grounds  that  finely  divided  powder  particles 
would  vibrate  as  if  each  was  a  separate  molecule,  thus 
giving  rise  to  an  extra  specific  heat  term  in  addition 
to  the  Montroll  correction.  This  contribution  was  ex¬ 
pected  to  approach  a  constant  value  at  low  temperatures, 
similar  to  Schroder's  superparamagnetic  contribution. 
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In  several  ensuing  experimental  programs  with  fine 
powders,  notably  the  work  of  Lien  and  Phillips  discussed 
earlier,  this  Jura  and  Pitzer  constant  value  term  was 
never  detected.  We  can  only  speculate  that  the  inter¬ 
particle  restoring  forces  postulated  by  them  were  too 
small  in  the  loose  powders  to  give  rise  to  the  required 
oscillations.  However,  in  the  case  of  a  metal  matrix 
containing  fine  cluster  structures,  such  as  those  considered 
by  Schroder,  such  restoring  forces  can  be  readily  imagined. 

Primarily  because  of  the  similarity  in  the  final  equation 
coming  from  Schroder's  classic  analysis  and  the  early 
theoretical  treatment  of  fine  powders  by  Jura  and  Pitzer, 
the  following  hypothesis  was  suggested: 

Hypothesis:  The  low  temperature  specific  heat 

of  certain  metal  alloys  containing 
small  ferromagnetic  cluster  structures 
can  be  treated  similar  to  a  matrix  of 
finely  divided  powder  and  thus  represent¬ 
ed  by  a  Montroll-type  equation  with  the 
addition  of  a  constant  term  representing 
the  contribution  of  cluster  oscillations. 
Consequently,  it  is  hypothesized  that  a  Montroll  equation 
of  the  type , 

C?ef-v<*T%(pe)>T  ('to) 

can  represent  these  metal  cluster  structures. 

A  recent  series  of  low  temperature  measurements  on 

several  metal  alloy  materials  which  evidence  such  a 

structure  was  recently  reported  by  Collings  and  Jelinek, 

22 

et  al  .  Following  Schroder's  approach,  they  fitted  the 
data  to  an  algorithm  similar  to  Eq.  39  •  The  value 
they  found  for  the  ^T^term  is  several  orders  of  magnitude 
smaller  than  the  remaining  two  terms.  Consequently,  we 
can  question  the  statistical  validity  of  its  determination 
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since  the  and  "B  terms  are  the  predominant  contributors 
to  specific  heat  at  the  lower  temperatures. 

Recognizing  that  the  term  comparatively  small 

in  this  temperature  region,  we  attempted  to  fit  the  Collings 
and  Jelinek  data  set  to  a  Montroll-type  equation  with  an 
additive  Jura  and  Pitzer  constant  term,  with  the  results 
(for  310S  stainless  steel): 


o 
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compared  with  the  authors'  statistical  fit  for  the  same 
material , 
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Our  statistical  fit  of  Eq.  4l  resulted  in  a  high 

correlation  coefficient.  Realistically,  however,  the 
.  •  2 

determination  of  the  T  term  is  open  to  the  same  criticism 
as  the  Collings  and  Jelinek  determination  of  the  T-^ 
coefficient  since  it  is  an  order  of  magnitude  smaller 
than  the  other  terms. 

Since  no  other  "raw"  data  has  been  located  in  the 
literature,  it  is  not  possible  to  further  test  this 
hypothesis  and  resolve  the  question  as  to  whether  the 
surface  enhancement  effect  of  interest  in  this  study  is 
automatically  associated  with  metal  alloy  cluster  structure 
materials.  The  production  and  handling  of  materials  which 
evidence  the  Montroll-type  enhanced  specific  heat  may 
obviously  be  considerably  simplified  if  cluster-forming 
alloys  which  meet  the  requirements  for  large  surface 
specific  heat  defined  by  Eq.  22  are  available. 

In  order  to  predict  whether  or  not  such  characteristics 
may  exist,  it  is  necessary  to  determine  the  effective  Debye 
temperature  for  alloys  which  form  this  type  of  cluster 
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structure.  Debye  temperature  tabulations  for  a  few  of 
such  alloys  have  been  located.  However,  since  a  broad 
range  of  transition  element  metals  can  form  alloys  of  the 
superparamagnetic  cluster  structure  type,  an  analytical 
method  must  be  utilized  to  predict  the  Debye  value  for  the 
most  likely  candidates.  The  development  of  this  method¬ 
ology  was  not  pursued  during  this  phase  of  the  study 
program. 


2 . 7  The  Composite  Specific  Heat  Model 

Since  all  Debye -type  specific  heat  models,  including 
Montroll's,  undercount  the  allowable  vibration  modes  of 
substances  exhibiting  a  large  surface-to-volume  ratio, 
the  computed  value  for  specific  heat  is  expected  to  be 
smaller  than  actual  measurement  results.  Reinforcing  this 
conclusion,  the  error  in  the  T2  (surface  contribution) 
term  of  the  Montroll  equation  has  been  found  to  range  in 
both  dielectric  and  metal  powders  from  a  factor  of  four 
to  nearly  an  order  of  magnitude  too  small,  based  on  the 
few  measurement  programs  reported  to  date. 

The  most  reliable  prediction  of  composite  specific 
heat,  then,  can  probably  be  obtained  by  an  empirical  modi¬ 
fication  to  Eq.  40,  The  suggested  form  iai 


C»  ST1  ♦‘■ton-*  ♦  •+ 15  .  CO) 

The  coefficients  in  this  equation  can  be  computed  as 
follows  I 


Cftu/|V\OW*/VH 


(20,  ref) 


0< 


(24,  ref) 
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(44) 


'B*  O(^)  CAk/houe/k7'  (45) 

The  notation  in  Eq.  45  indicates  that  the 
coefficient  ""5  (associated  with  superparamagnetic  metal 
alloys  only)  has  a  numerical  value  on  the  order  of  the 
linear  coefficient,  gamma.  This  approximation  is 
based  on  limited  test  data  stemming  from  the  Collings  and 
Jelinek  work,  and  may  not  prove  to  be  generally  valid. 

There  apparently  has  been  no  successful  theoretical  develop¬ 
ment  of  an  equation  to  predict  Das  yet,  since  no  work 
along  these  lines  was  discovered  in  the  literature.  It 
was  noted  earlier  that  cluster  structure  materials  of 
this  kind  were  not  considered  in  depth  in  this  phase  of 
the  study;  consequently,  there  was  no  attempt  as  part  of 
this  program  to  model  an  equation  for  this  term. 

2.7.1  Comparison  with  Reported  Measurements  -  In  order  to  get 
some  indication  of  the  likely  accuracy  of  Eq.  43  , 
its  specific  heat  predictions  for  dielectric  materials 
were  compared  to  the  Lien  and  Phillips  data  on  K {O  and 
the  Barkman,  et  al  data  on  NxvCfi. .  Results  are  shown  in 
Figures  3  and  4  respectively.  Similarly,  the  com¬ 
puted  predictions  are  compared  to  the  palladium  data  of 
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T2  (°K2) 

figure  3.  Magnesium  Oxide  Small  Particle  Cv 
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Comsa,  et  al ,  in  Figure  5  and  the  Novotny  and  Meincke 
data  for  lead  and  indium  in  Figures  6  and  7* 

The  Montroll  theory  and  the  empirical  equation  devel¬ 
oped  earlier  in  this  section  for  surface  enhanced  specific 
heat  can  be  correlated  reasonably  well  with  these  few 
reported  measurement  sets  except  for  the  1976  data  of 
Novotny  and  Meincke  on  lead  and  indium  particles.  The 
correlation  in  this  case  is  poor  if  we  assume  that  the 
Eq.  43  specific  surface  value  (which  wasn't  measured 
by  Novotny)  is  the  same  as  it  would  be  for  spherical 
particles.  A  failure  of  the  theory  for  these  materials 
would  be  particularly  serious,  since  they  both  are  low 
Debye  temperature  metals,  and  prime  candidates  for  the 
regenerator  improvement  test  program. 

The  Novotny  measurement  results  were  derived  under 
unusual  test  conditions,  which  may  have  invalidated  the 
conditions  necessary  for  Eq.  43  to  hold.  To  assess 
the  probable  influence  of  these  test  conditions  on  the 
reported  measurements,  it  is  necessary  first  to  examine 
the  more  recent  data  set  of  Comsa,  et  al ,  performed  on 
palladium  ultrafines.  This  palladium  data  fits  the 
predictions  derived  from  the  Montroll  theory  exceedingly 
well . 

Comsa  and  his  co-workers  prepared  the  ultrafines  by 
the  inert  gas  evaporation  technique,  allowing  the  palla¬ 
dium  vapor  to  condense  in  a  helium  environment  to  minimize 
particle  size.  Granqvist  and  Buhrman  ^  recently  showed 
that  sub-200  ultrafines  produced  in  this  manner  are  spher 
ical  (as  opposed  to  taking  the  characteristic  crystalline 
shape),  and  have  a  narrow  spread  in  their  size  distribution 
characterized  by  a  log  normal  distribution  function  (LNDF). 

Since  Eq.  43  strictly  applies  to  a  monodisperse , 
or  single  size,  distribution  of  particles,  it  is  apparent 
that  any  measurable  quantity  of  ultrafines  larger  than  the 
nominal  size  will  invalidate  the  analytical  prediction  of 
specific  heat.  The  larger  particles  will  contribute  less 
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of  a  specific  heat  enhancement  than  predicted,  and  take 
up  a  disproportionate  share  of  the  sample  mass,  causing 
the  analysis  to  err  both  in  the  increment-per-particle 
to  be  expected  as  well  as  in  the  estimate  of  total  number 
of  particles  which  give  this  contribut ion . 

The  Comsa  measurement  program  included  electron  micro¬ 
graph  studies  of  the  palladium  ultrafines  which  seems  to 
confirm  the  narrow  size  distribution  predicted  by  Granqvist. 
Based  on  this  palladium  data,  we  conclude  that  the  LNDF 
distribution  is  acceptably  narrow  for  purposes  of  achiev¬ 
ing  the  full  benefit  predicted  by  the  Montroll  analysis. 

In  contrast  to  the  Comsa  program,  Novotny  prepared 
his  samples  by  impregnating  a  porous  glass  with  metal 
vapors.  His  program  was  carefully  designed  to  determine 
the  diameter  of  the  capillary  pores,  but  not  their  depth; 
his  size  distribution  of  the  pore  diameters  appears  to  be 
relatively  narrow. 

The  resulting  metal  particle  shape,  similar  to  the 
glass  pore  shape,  apparently  tends  toward  long  cylinders 
instead  of  discrete  separated  spheres,  and  this  could 
explain  his  conclusion,  contradictory  to  Comsa,  that  the 
crystal  habit  is  evident  in  the  particles. 

Novotny  finds  evidence  of  a  significant  surface  en¬ 
hancement  effect,  but  the  measured  magnitude  does  not 
follow  the  prediction  of  Eq.  43  as  applied  to  spherical 
particles  of  the  pore  diameter.  The  discrepancy  is 
expected  due  to  the  fact  that  the  long  cylinders  of  metal 
in  the  glass  capillary  pores  have  a  different  specific 
surface  than  would  spheres  of  the  same  diameter.  Eq.  43 
fits  the  data  of  Figure  6  reasonably  well  if  we  assume 
that  the  22$\  diameter  metal  cylinders  approximate  spheres 
of  about  130  -  150&  diameter.  For  the  22A,  indium  particles, 
Novotny's  results  can  be  fit  reasonably  well  under  the 
assumption  that  the  equivalent  particle  size  is  in  the 

O 

80  -  100  A  range,  implying  that  the  cylinders  of  indium 
in  the  pores  are  somewhat  shorter  than  those  of  the  lead. 
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If  Novotny  had  measured  the  specific  surface  of  his 
samples,  thus  taking  into  account  their  deviation  from 
a  simple  spherical  shape,  then  Eq.  24  could  have 
been  used  directly  without  making  an  assumption  of 
spherical  structure.  We  can  only  speculate  that  the 
resulting  correlation  to  his  measurements  might  have  been 
significantly  better. 

In  summary  then,  the  Novotny  results  do  indeed  detect 
a  significant  surface  enhancement  in  both  the  indium  and 
lead  ultrafine  cylinders  produced  by  the  porous  glass 
test  technique.  Some  combination  of  shape  factor  and/or 
non-homogeneity  in  the  pore  size  distribution  probably 
explains  the  fact  that  specific  heat  measurement  of  these 
cylinders  deviates  from  the  results  predicted  for  simple 
spheres.  In  contrast  to  Novotny's  work,  the  recently 
reported  palladium  data  shows  convincingly  that  spherical 
ultrafines  can  be  produced  (by  the  inert  gas  evaporation 
technique)  and  that  the  resulting  size  distribution  is 
narrow  enough  that  the  Montroll  approach  reasonably 
accurately  predicts  the  specific  heat  enhancement. 

In  order  to  verify  that  the  "long  cylinder"  theory 
adequately  explains  the  Novotny  results,  another  set  of 
measurements  using  spherical  ultrafines  of  lead  and 
indium  will  have  to  be  performed. 

2 . 8  Candidate  Materials  for  Surface  Specific  Heat  Enhancements 

According  to  Eq.  25  >  the  most  important  selection 
criterion  for  materials  likely  to  show  a  large  surface 
enhanced  specific  heat  effect  in  the  cryogenic  region  is 
the  bulk  Debye  temperature.  Of  lesser  importance  are  the 
bulk  density  and  molecular  weight  of  the  candidate  sub¬ 
stances. 

Since  the  Debye  temperature  actually  varies  slightly 
with  temperature  in  the  region  where  surface  effects 
become  important,  the  low  temperature  value ,  Qo  ,  has 
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been  used  in  the  following  tabular  summary,  where  available, 
instead  of  the  conventional  value  of  ©0.  The  difference 
in  the  two  is  that  is  evaluated  close  to  absolute  zero 
while  the  latter  is  more  an  average  value  evaluated  in 
the  vicinity  of  3/2  . 


To  illustrate  the 

differences  which 

can  be  found, 

2  4 

Table  1  lists  some  typical  comparisons 

©o 

©o 

Bulk  Material 

( average) 

(low  temperature) 

Aluminium  (Al) 

385 

426 

Lead  (Pb) 

85 

88 

Indium  (In) 

140 

109 

Palladium  (Pd) 

275 

262 

Rubidium  (Rb) 

60 

55 

Mercury  (Hg) 

100 

72 

Table  1.  Variation  in  %  and  ©o  for 
Various  Substances 

Table  2  presents  pertinent  data  for  nineteen 
candidate  materials  which  have  been  found  to  exhibit 
favorable  values  of  Debye  temperature. 

Eight  of  these  materials  appear  promising  enough 
at  first  inspection  that  plots  of  molar  specific  heat  for 
each  are  presented  in  Figures  8  through  15,  and  volumetric 
specific  heat  in  Figures  16  through  23.  These  candidates 
include,  in  order,  bismuth,  mercury,  indium,  lanthanum, 
lead,  selenium,  thallium,  and  chromium  chloride.  Each 
figure  illustrates  calculated  specific  heat  values  in 
the  5  -  20°K  range  as  determined  from  Eq.  43  for  the  bulk 
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Figure  8.  Molar  Specific  Heat  of  Bismuth  Ultrafines 
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Figure  9.  Molar  Specific  Heat  of  Mercury  Ultrafines 
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Figure  10.  Molar  Specific  Heat  of  Indium  Ultrafines 
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Figure  11.  Molar  Specific  Heat  of  Lanthamuro  Ultrafinea 
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Figure  12. 
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Molar  Specific  Heat  of  Lead  Ultrafines 
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Figure  13.  Molar  Specific  Heat  of  Selenium  Ultrafines 
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MOLAR  SPECIFIC  HEAT 


Figure  14,  Molar  Specific 
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Figure  15.  Molar  Specific  Heat  of  Chromium  Chloride  Ultrafines 


f  o  o 

material,  100  A*,  ,  50  A  and  22  A,  ultrafines,  plotted 

against  bulk  lead  as  the  baseline  reference. 

Upon  inspection  of  these  curves  CeCW  can  be  eliminated 
from  further  consideration  because  of  its  low  density,  and 
mercury,  lead,  and  selenium  remain  as  the  primary  candidates 
for  the  regenerator  improvement  experiment  to  be  conducted 
in  later  phases  of  this  program. 


2.8.1  Interpretation  of  Specific  Heat  Enhancement  Data.  It  was 
noted  in  Eq.  15  that  the  ^"^approximation  for  bulk 
material  specific  heat  is  strictly  valid  only  at  temper¬ 
atures  much  lower  than  the  characteristic  Debye  temperature 
of  the  solid.  At  temperatures  greater  than  about  ©tM 
the  approximation  shouldn't  be  used  for  high  precision 
work,  and  at  temperatures  warmer  than  about  3>/a  the 
complete  form  of  the  Debye  equation: 


C,,~"  *’J0  (e*-  4 


©D 


:  *" 


(46) 


yields  lower  specific  heat  values  than  does  the 
approximation  (Figure  24),  becoming  increasingly 
significant  as  T  increases. 

All  measurements  done  on  fine  powders  to  date  have 
been  at  temperatures  where  the  simplified ^T\pproxi- 
mation  for  the  bulk  material  applies.  The  data  of 
Figures  3  and  4  on  t^O  and  N*Cl  shows  clearly  that  a 
significant  surface  enhancement  in  specific  heat  occurs 
at  these  fractional  Debye  temperature  regions. 
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Figure  22.  Volumetric  Specific  Heat  of  Thallium  Ultrafines 
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Figure  23.  Volumetric  Specific  Heat  of  Chromium  Chloride  Ultrafines 
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The  recent  measurements  on  palladium  ultrafines 
illustrate  that  the  surface  enhancement  effect  does 
occur  similarly  with  metals,  and  further  shows  that 
the  empirical  form  of  the  algorithm  (Eq.  43  )  produces  a 
conservative  fit  in  this  fractional  region. 

The  warm  end  of  these  reported  palladium  measure¬ 
ments  approached  a  temperature  value  of  .  There 

is  no  indication  at  that  temperature  that  the  enhanced 
specific  heat  is  beginning  to  deviate  from  the  value 
predicted  by  Eq.  43. 

The  Montroll-type  theory  isn't  developed  far  enough 
to  indicate  the  expected  warm-end  temperature  limits  on 
the  enhancement,  and  no  reported  measurements  have  in¬ 
vestigated  this  effect  beyond  the  range  of  the  palladium 
data. 

In  order  to  apply  this  surface  specific  heat  theory 
to  cryogenic  refrigerator  regenerator  design  it  is  necessary 
to  predict  the  enhancement  effect  at  temperatures  as  warm 
as  about  .  Because  Montroll's  equation  was  derived 

under  the  limit  T— «*0,  we  expect  thee*T**  correction  to 
be  in  error  (similar  to  Debye ' s  term)  as  the  value 
of  ^>/r  becomes  smaller.  Since  no  test  data  exists  we 
must  either  assume  that  Eq.  43  is  valid  all  the  way 
out  to  the  warm  end  of  our  interest,  or  else  artificially 
assume  some  degradation  in  the  computed  value  between 

©b/g 

(where  palladium  data  is  available)  and^/"&  .  So 
as  to  attempt  to  keep  performance  projections  conservative, 
we  shall  assume  that  the  full  calculated  value  of  the 
enhancement  increment  is  applicable  as  warm  as  &p/%  --as 
indicated  by  the  palladium  data  --  and  decreases  percent¬ 
agewise  to  zero  (in  a  manner  linear  with  temperature)  at 
a  value  of  .  Figure  23  illustrates  this  hypothesized 

fall-off  in  the  c<Tl  term  for  our  baseline  material  (lead). 
Specific  heat  values  from  this  figure  will  be  used  in 
Section  III  when  enhanced  refrigerator  performance  predictions 
are  compared  with  performance  derived  from  conventional 
technology  using  lead  spheres  in  the  regenerator. 


Figure  24 .  Comparison  of  Specific  Heat  Calculations  for  Lead 
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IDEALIZED  REFRIGERATOR  PERFORMANCE 


3.1  Introduction 

The  mechanisms  used  to  produce  refrigeration  at 
cryogenic  temperatures  have  historically  been  inefficient 
devices,  plagued  with  reliability  problems  and  high  input 
power  requirements.  In  the  early  1960s  several  research 
groups  began  showing  promising  results  in  applying  the 
refrigeration  version  of  the  Stirling  heat  engine  to  the 
production  of  low  temperatures.  Used  in  this  manner, 
the  process  is  sometimes  called  the  Kirk  cycle,  after  an 
early  researcher  who  reported  its  utilization  to  produce 
refrigeration.  Following  the  more  usual  nomenclature, 
however,  it  will  be  referred  to  as  the 
Stirling  cycle  in  this  report,  and  is  shown  graphically 
in  Figure  26. 

The  second  highly  significant  development  came  in 
that  same  time  frame  when  Gifford  and  McMahon  published 
a  report  describing  a  new  thermodynamic  cycle  which  now 
bears  their  names  J .  There  has  previously  been  a  debate 
whether  their  disclosure  actually  represented  a  new  thermo¬ 
dynamic  concept  or  was  merely  a  novel  implementation  of 
the  Ericsson  cycle  shown  in  Figure  27.  Regardless  of 
how  we  argue  the  resolution  of  that  question,  the 
implementation  of  their  conoept  resulted  in  high  reliability 
and  high  performance  machines  when  compared  with  prior 
technology. 

The  Stirling  cycle  and  the  Gifford-MoMahon  cycle 
refrigerators,  as  well  as  devices  based  on  several  other 
similar  cycles  not  highlighted  here,  have  in  common  the 
fact  that  they  are  designed  around  a  periodic  flow  heat 
exchanger  called  a  regenerator.  The  regenerator  separates 
the  hot  working  volume  of  these  machines  from  the  cold 
working  volume,  and  is  the  keystone  of  machine  performance. 
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Figure  26  .  Stirling  Cycle  Diagram 


An  ideal  regenerator  would  have  characteristics  of  zero 
pressure  drop  for  gas  flow  between  the  hot  and  cold 
regions,  infinite  specific  heat  so  as  to  cause  a  perfect 
temperature  transition  (between  hot  and  cold)  for  the  gas 
flowing  through  it,  zero  heat  conductivity  along  the  axis 
of  the  heat  exchange  matrix,  and  an  infinitesimally 
small  void  volume  in  the  gas  flow  paths.  Obviously, 
these  ideal  requirements,  even  to  the  extent  they  are 
not  mutually  exclusive,  cannot  be  satisfied.  Consequently, 
the  performance  of  the  overall  refrigerator  is  directly 
affected  by  the  deviation  of  the  regenerator  from  the 
ideal . 

Because  the  Stirling  cycle  machine  is  subjected  to 
relatively  severe  mechanical  loading  on  dynamic  components 
such  as  bearings  and  seals,  a  closely  related  thermo¬ 
dynamic  process  which  minimizes  these  loads,  known  as 
the  Vuilleumier  (VM)  cycle,  was  adopted  for  these  refrig¬ 
erators  in  the  late  1960s.  VM-based  machines  have 
incrementally  improved  the  reliability  of  Stirling- type 
devices,  and  are  leading  candidates  for  long  life  space 
mission  applications. 

Some  writers  describe  the  VM  cycle  as  merely  a 
Stirling  cycle  modified  so  that  the  high  pressure  is 
produced  by  a  thermal  (rather  than  a  mechanical)  compressor. 
It  is  probably  more  accurate,  however,  to  describe  it  as 
a  variant  of  the  Ericsson  cycle,  since  mass  flow  out  of 
the  cold  volume  occurs  during  the  expansion  part  of  the 
cyele2^.  The  thermodynamic  process  which  is  physically 
implemented  in  the  VM  refrigerators  cannot  conveniently 
be  described  on  a  temperature -entropy  plot,  because 
dynamic  flow  conditions  are  encountered  throughout  the 
cycle.  Flow  process  thermodynamics  are  usually  employed 
to  analyze  the  performance  of  these  machines, 

The  first  law  of  thermodynamics,  for  steady  flow 
processes  through  a  control  volume,  is  written  as: 
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when  kinetic  and  potential  energy  terms  are 
neglected . 


QvOeVi 


is  net  heat  transfer  rate  to  or  from  the 


Wviefi 

H.  o  i 1 

jC/wHji  i: 


net  work  rate  done  by  the  system 
the  enthalpy  sum  for  all  flow  outlets 
the  enthalpy  sum  for  all  flow  inlets. 


system 


It  is  assumed  that  over  the  period  of  the  process  being 
analyzed  there  is  no  net  storage  of  mass  or  energy  in 
the  control  volume  considered. 

This  approach  is  adapted  to  the  analysis  of  the 
periodic  low  associated  with  the  expansion  volumes  of 
the  regenerator-based  refrigerators,  and  the  equation  for 
gross  refrigeration  capacity  per  cycle  is  usually  expressed 
as : 

""  (48) 

where 


<sw  indicates  the  heat  absorbed  in  the  actual 
'  expansion  process  (in  contrast  to  the  Qj  ^ 
idealized  quantity  from  Figure  26a) 
is  the  work  done  by  the  expanding  gas,  and 

is  the  enthalpy  change  of  the  gas  during 
the  expansion  process. 


Helium  is  almost  always  employed  as  the  working  gas 
in  refrigerators  of  this  type,  since  it  liquefies  only 
at  pressures  and  temperatures  lower  than  those  normally 
achieved  in  the  cycle.  For  purposes  of  computing  gross 
refrigeration,  the  assumption  is  often  made  that  helium 
behaves  as  an  ideal  gas,  following  the  relationship: 


Tv  ■=  ™"£rr 


(49) 


where  /wv.  is  the  mass  of  gas  involved  in  the 
process,  and  is  the  gas  constant. 

For  a  closed  system,  the  first  law  of  thermodynamics  is, 
in  analogy  to  Eq.  2  : 

cal  U-  avJ  (50) 

where  the  symbols  indicate  that  d  l)  is  an  exact 
differential  whose  value  depends  only  on  the  state  of 
the  thermodynamic  system,  while  SQ  and  are  inexact 
differentials  and  represent  infinitesimal  quantities 
whose  values  depend  on  the  particular  process  followed 
in  taking  the  working  substance  from  state  to  state. 
Since 


S  MO- 'PdV 


we  can  write 

q|U  ” 


(5D 


(52) 


and  for  a  constant  volume  process,  we  showed  in  Eq.  5 
that 


(53) 


since  the  internal  energy  of  an  ideal  gas  is  a  function 
of  temperature  only,  i.e.: 

0=U(t).  (5*0 
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Since  we  define  enthalpy  as 


tt  =  U+fV  , 

we  can  write  for  an  ideal  gas 

hi  =  L)  ^  yw\^*T  ^ 


(55) 


(56) 


and  since  all  terms  on  the  right  side  of  the  equation 
are  solely  functions  of  temperature,  it  follows  that  H  , 
likewise,  is  only  a  function  of  temperature. 

Consequently,  if  the  control  volume  selected  for 
the  regenerator  flow  analysis  includes  both  the  regener¬ 
ator  and  the  expansion  volume ,  and  it  is  postulated 
that  inlet  and  exit  gas  are  at  the  same  temperature 
(i.e.,  a  100$  efficient  regenerator),  then  the  equation 
for  gross  refrigeration  can  be  rewritten  as: 


Q'  =  \  TW  <5?> 

where  the  prime  ( i )  emphasizes  that  the  ideal  gas 
assumption  was  employed  to  achieve  this  simplification, 
and  the  symbol  ^.indicates  that  the  work  term  is  in¬ 
tegrated  over  the  expansion  cycle. 

In  summary  then,  for  VM  cycle  refrigerators  the 
estimate  for  gross  refrigeration  production  uses  either 
Eq.  48  or  the  ideal-gas  simplification,  Eq.  57  ■ 

As  a  matter  of  fact,  since  all  the  cycles  of  interest  as- 
implemented  deviate  from  the  ideal  'P—V  diagrams,  these 
two  equations  can  be  applied  to  any  of  these  machines  to 
estimate  the  actual  output. 

The  first  successful  attempt  to  analyze  the  performance 
of  non-idealized ,  as -implemented  Stirling  machines  was 
reported  by  Schmidt  in  the  l860s2^.  His  analysis  recog¬ 
nized  that  realizable  machines  would  employ  sinusoidal 
drives  for  the  reciprocating  elements  rather  than  the 
"square  wave"  drives  implied  in  the  piston  movements  of 
Figure  26. 
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His  theory  requires,  among  other  simplifications, 
the  assumptions  of  isothermal  compression  and  expansion, 
and  further  assumes  the  regenerator  operates  with  lOOfo 
efficiency.  In  spite  of  the  fact  that  his  model  is  still 
highly  idealized,  it  represented  a  striking  improvement 
over  the  cycle  analysis  equations  based  on  the  pure  "P— V 
diagram  for  the  refrigerator. 

This  theory  has  been  further  refined  in  recent  years. 
The  Philips  Laboratories,  Eindhaven,  Netherlands  have 
been  heavily  engaged  in  research  and  development  of 
Stirling  devices  for  several  decades.  Early  emphasis 
was  on  the  Stirling  engine  as  a  prime  mover;  however, 

I 

as  early  as  1945  Rinia  and  duPre  of  that  organization 
were  involved  in  preliminary  work  on  refrigerators. 

More  recently,  computer  assisted  design  analyses  of  VM 
refrigerators  have  been  published.  Correlation  between 
as-built  performance  and  the  computer-predicted  capability 
is  reportedly  good  in  at  least  one  case  of  a  high  capacity 

o  Q 

VM  machine  developed  for  the  Air  Force^  . 


3-2 


Stirling  Cycle  Analysis 

For  ease  in  understanding  the  concepts  to  be  presented 
in  this  report  --  specifically,  the  concept  of  improving 
overall  refrigerator  performance  by  a  technique  to  increase 
the  specific  heat  of  regenerator  matrix  materials  --  it 
is  particularly  convenient  to  avoid  both  the  Schmidt 
analysis  and  the  flow  process  analysis  and  to  work  with 
a  true  cycle  which  can  be  described  on  a"PV  state  diagram. 
Consequently,  the  Stirling  cycle  of  Figure  26a  will  be 
used  hereafter  as  our  model  for  discussing  this  effect. 

The  idealized  cycle  achieves  the  same  theoretical 
coefficient  of  performance  (COP)  as  the  theoretical 
Carnot  cycle,  whose  limiting  performance  forms  an  alternate 
definition  of  the  second  law  of  thermodynamics.  COP  is 
defined  as: 


(58) 
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Where  the  subscript  means  idealized  (or  Carnot)  COP, 

means  the  gross  heat  absorbed  in  going  from 

State  3  to  4  of  Figure  26,  and  V\/v,Q+  is  the  difference 

ne  u 

between  the  theoretical  work  of  compression  at  the  high¬ 
er  temperature  and  the  expansion  work  available.  For  the 
idealized  cycles 

(59) 

T*h 

Before  proceeding  to  consider  the  effect  of  regenerator 
efficiency  on  machine  performance  as  well  as  the  effect 
of  adding  multiple  stages,  it  is  useful  to  list  the  heat 
and  work  balance  equations  for  the  cycle.  For  the 
isothermal  compression  1-2  we  can  apply  the  first  law 
from  Eq.  2  as: 


Q»-z~  ^1-z.  (60) 

and  our  sign  convention  requires  that  Qbe  numerically 
positive  for  heat  absorbed  by  the  system,  and  ^  be 
numerically  positive  for  work  done  by  the  system. 

Treating  the  helium  as  an  ideal  gas 

H-i '  a* 

thus 

qux--  rth  !L % 

For  the  isochoric  process  2-3.  the  first  law  becomes, 


(60a) 


(60b) 


S  ^Z-3 

and  consequently 


(60c) 
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4 


( 60d) 


and,  of  course 


With  similar  arguments, 

~  w'UT’c  JU  *<4* 

and 

<V, 3  ^cv  CYTc) 


while 


oo 


4-i 


-  o 


We  thus  find  that, 


(T.-T„  ^  . 


( 60e ) 


(60f ) 


(60g) 


( 60h) 


(61) 


(6la) 


3*2.1  Simple  Regenerator  Model.  In  the  ideal  cycle  of 

Figure  26,  the  heat  transfer  in  both  Steps  2-3  and 
4-1  is  accomplished  in  the  high  efficiency  periodic 
flow  heat  exchanger  referred  to  earlier  as  a  regenerator. 

In  the  figure  the  entire  mass  of  helium  contained  in  the 
cycle  is  required  to  be  transformed  exactly  between  the 
two  temperature  levels  as  it  passes  through  the  regenerator. 
In  real  machines  it  is  not  possible  to  achieve  such  per¬ 
fect  heat  transfer,  and  consequently  the  gas  exiting  from 
the  cold  end  of  the  regenerator  will  be  warmer  thanl^  by 
some  &T. 
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We  can  define  regenerator  efficiency  (£)  in  terms 
of  the  heat  actually  transferred  during  the  flow  proces 

_  &C.4- 

6*_2lL_  (62^ 

Q* 

where  the  superscript  act  refers  to  the  heat 
actually  transferred  during  the  process  and  id 
refers  to  the  idealized  amount  of  heat  which 
should  have  been  transferred  to  properly  "fit" 
the  cycle  diagram.  We  know  from  the  earlier 
derivation  that: 

(T£-TH^  .  (63) 

if  we  define  the  regenerator  heat  loss  as 

(64) 


we  can  show 


Qp*  6’^)mCv 


It  is  also  necessary  to  recognize  that  in  a  real 
machine  the  gross  refrigeration  produced  in  the  process 
3-4  is  "used  up"  by  the  sum  of  (A)  the  external  load 
being  cooled  (  )f  (B)  the  regenerator  heat  loss 

(  3*  ),  and  (C)  various  other  stage  losses  associated 

with  the  design  of  the  particular  machine  of  interest 
(  Ovu.  ).  This  can  be  expressed  as: 


(66) 


We  saw  earlier  that, 


This  can  be  rearranged  as  _ 

where 


Then 

Q«  _ 

Q>4 


0-6)  C^/ro  -0 

(jr‘*  JWMO 


(67) 

(68) 

(69) 

(70) 


Equation  70  can  be  solved  to  determine  that  minimum 
efficiency  which  allows  all  the  gross  refrigeration  to 
be  consumed  by  the  regenerator  heat  loss  (this  approach 
can  be  modified  to  include  Q>  L>- 
Thus : 

-  l “ 

A  graph  of  Eq.  71 
limiting  case  there  is  no  net  refrigeration  left  to  cool 
an  external  heat  load. 


(r*Ai,  ~  O  <7 

is  included  as  Figure  28.  In  this 


FOR  HELIUM,  ^  s  1.4,  TH=  324#K,  V=  1.4 


Figure  28.  Minimum  Required  Regenerator  Efficiency 
re.  Tc  (Cold  Stage) 


(72) 


Now,  the  theoretical,  or  design,  efficiency  €,©  of  a 
periodic  flow  regenerator  can  be  expressed  as: 


X. 


where  A  o  is  the  theoretical  number  of  heat 
transfer  units  associated  with  the  regenerator 
design . 

The  value  for  the  number  of  design  heat  transfer  units 
is  computed  as: 


(73) 


where  frf  is  the  heat  transfer  surface  area 

is  the  overall  coefficient  of  heat 
transfer,  approximated  for  our  purposes 
by  the  convection  heat  transfer 
coefficient ,  M. , 

and  C*  is  the  gas  flow  stream  capacity  rate 

*  A 

=  /*v\  Cp. 

We  can  thus  write  the  approximation. 


\  J*  frill  (71*) 

6  “ 

Since  the  convection  coef f icient  ,  is  not  always 

accurately  computed  from  basic  heat  transfer  principles, 
it  is  useful  to  determine  it  indirectly  for  a  particular 
heat  exchanger  configuration  from  experimental  evaluation 
of  the  Stanton  number,  >  a  dimensionless  heat  transfer 

modulus  defined  as: 
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&_ 

Kt=  ”5' 

where ^  is  the  regenerator  flow-stream  mass 
velocity, 

<761 

with  Kc  being  the  minimum  free -flow  area  of 
the  heat  exchanger. 

If  we  neglect  entrance,  exit,  and  flow  acceleration 
effects  in  the  regenerator,  Fanning's  equation  which 
describes  the  flow  pressure  drop  due  to  core  friction 
can  be  written  as: 


tv  ■£  6J_L 

T  =  *-k 

(77) 

where 

&P  is  core  pressure  drop 

and 

L.  is  length  of  the  flow  path 

Jf-  is  the  Fanning  friction  factor 
^  is  average  gas  density 

is  the  hydraulic  radius  of  the 

defined  as 

flow  path 

rt 

(78) 

If  we  define  the  heat  transfer  "J"  factor  as, 


where  Nfj^is  the  Prandtl  number: 


(79'1 


2 


(80) 


78 


then  Fanning's  equation  can  be  rewritten  to  relate  the 
pressure  drop  and  heat  transfer  characteristics  of  the 
regenerator  in  the  form: 

AP  <a*£ 

X’  *  J  r  * 


(81) 


Consequently,  the  required  length  of  the  regenerator 
heat  transfer  matrix  can  be  determined  by  combining  the 
two  forms  of  Fanning's  equation. 


L 


pa 


(82) 


3.2.2  Multiple  Stage  Refrigeration.  Up  to  this  point  we  have 
been  considering  a  simple  "single  stage"  refrigerator 
wherein  a  single  regenerator  was  used  to  span  the  temper¬ 
ature  range  between  and  tc. 

It  has  been  found  experimentally,  however,  that 
machines  capable  of  achieving  temperatures  lower  than 
50  -  ?0°K  achieved  their  maximum  efficiencies  if  several 
stages  of  cooling  were  employed.  This  staged  cooling 
allows  heat  leaks  to  be  intercepted  at  temperatures 
warmer  than  the  desired  design  temperature.  In  addition, 
the  stage  concept  allows  extra  refrigeration  to  be  made 
available  (for  example,  to  cool  radiation  shields  at 
these  intermediate  stage  temperatures)  with  a  relatively 
modest  penalty  in  work- input,  compared  with  providing 
that  same  quantity  of  refrigeration  at  the  lowest  design 
temperature . 

In  order  to  understand  the  performance  of  a  typical 
three -stage  refrigerator  a  model  can  be  developed  where¬ 
in  the  gas  expanded  at  each  stage  undergoes  a  series 
of  processes  which  each  can  be  described  on  a  P-V  diagram 
as  illustrated  in  Figure  29.  Each  volume  of  gas  works 


Figure  29.  P-V  Diagram  Representation  for  Multiple  Stage  Refrigerator 
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AC{1)  :  FLOW  CROSS  SECTION  FOR  STAGE  1 
AC(2)  :  FLOW  CROSS  SECTION  FOR  STAGE  2 
Aq(3)  :  FLOW  CROSS  SECTION  FOR  STAGE  3 


between  temperature  limits  established  by  the  cold  work¬ 
ing  temperature  of  that  stage  and  some  established 
typically  by  the  crankcase  ambient  temperature  of  the 
machine . 

The  flow  path  for  each  stage  is  conveniently  con¬ 
ceptualized  as  shown  in  Figure  30.  Neglecting  mixing 
(which  is  immaterial  as  long  as  we  have  "perfect"  re¬ 
generators),  the  gas  to  be  expanded  at  the  first  (warmest) 
refrigeration  stage  can  be  considered  to  flow  through  the 
length  of  the  first  stage  regenerator,  requiring  some  free¬ 
flow  area  't^e  overall  regenerator  frontal  area. 

Likewise  the  gas  to  be  expanded  in  Stages  2  (intermediate 
stage)  and  3  (coldest  stage)  will  require  the  amount  of 
flow  area  designed  as  ftcCx)  and  respectively.  We 

can  express  the  total  required  free  flow  area  of  this 
first  stage  regenerator  as: 

_co  ^ 

\  *  ^CCO  *  (83) 

Similarly,  the  total  free  flow  area  of  the  second  stage 
is, 

(84) 


(85) 

since  flows  earmarked  for  the  warmer  stages  undergo 
their  expansion  at  that  earlier  stage  and  then  flow  back 
to  the  ambient  temperature  crankcase. 

It  may  be  useful  to  think  of  this  multi-stage  process 
as  three  devices  operating  in  parallel  (not  in  series), 
each  one  independently  undergoing  the  processes  shown 
in  Figure  29 »  only  coincidentally  thermodynamically  coupled 


and  for  the  third  stage , 

2  *«t»)  j 
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via  the  common  mechanisms  of  the  refrigerator.  To 
further  belabor  this  point,  the  overall  regenerator 
temperature  gradient  required  for  proper  heat  transfer 
between  and  during  the  two  isochoric  processes 
of  the  cycle  is  provided  in  part  by  the  first  stage 
regenerator  (fromTTt  toTtl  ),  in  part  by  the  second 
stage  regenerator  (fromT^  toT^  ),  and  finally  by 
the  third  stage  regenerator  (from“^  to  T<.  )•  It 
is  only  in  this  sense  that  we  can  talk  about  the  stages 
being  arranged  in  "series"  in  this  particular  model. 

Equations  60  through  85  must  be  evaluated 
for  each  refrigeration  stage  during  the  process  of 
designing  a  multi-stage  machine.  For  purposes  of  this 
study  program,  however,  our  chief  interest  focuses  only 
on  the  variation  in  performance  of  the  coldest  stage  of 
such  a  device  caused  by  variations  in  the  specific  heat 
of  the  material  used  as  the  heat  storage  matrix  in  the 
regenerator . 

Recognizing  that  no  stage  of  a  multi-stage  refriger¬ 
ator  operates  as  an  ideal  heat  pump,  still  we  lose  very 
little  in  generality  by  assuming  that  whatever  the 
regenerator  efficiency  and  stage  losses  associated  with 
the  first  and  second  stages  of  a  particular  device  happen 
to  be,  they  do  not  change  with  a  variation  in  third  stage 
matrix  specific  heat,  and  therefore  treating  the  two 
warmer  stages  as  having  perfect  regenerators. 

Having  made  this  assumption,  the  work  input  necessary 
to  absorb  the  heat  loads  of  these  warmer  stages  can  be 
computed  based  on  Carnot  performance  then  merely  scaled 
up  in  accordance  with  the  relationship  a  manufacturer's 
typical  machine  has  to  the  theoretical  minimum  power 
requirements . 

3-2.3  Cold  Stage  Analysis.  When  we  focus  on  the  performance 

of  the  third  stage  regenerator,  we  find  that  its  actual 
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efficiency  (£)  differs  from  the  theoretical  design 
efficiency  ( ^0) ,  being  a  strong  function  of  the  ratio  of 
total  matrix  (thermal)  stream  capacity  (Cjj)  to  the  gas 
(thermal)  stream  capacity  (C^)>  The  matrix  stream 
capacity  is  defined  as  follows: 


where 

M 


I  ^  <**' 


(86) 


is  the  total  mass  of  heat  transfer  packing  trans- 
versed  by  the  gas  stream  during  the  cycle,  and  is 
equal  to  twice  the  regenerator  mass  (/VAg^) 
is  the  cycle  time  required  for  the  system  to  pass 
through  all  the  states  of  the  P-V  diagram,  and 
the  functional  notation  emphasizes  the  fact 

that  the  cold  stage  matrix  specific  heat  can  vary 
with  axial  location. 

We  can  also  define  blow  time  (^g)  as  the  time  required  for 
the  helium  gas  to  pass  through  the  regenerator  in  either 
direction,  and  approximate  this  time  by  the  relationship 


h:.  =  4b:, 


Thus 


<1  I— 


(87) 


(88) 


If  we  make  the  further  approximation  that  the  specific 
heat  of  helium  gas  is  constant  In  this  temperature  region, 
then  the  gas  stream  capacity  relationship  is 


s 


a 

-  0\!C 


(89) 


and  consequently, 

c* 


M 


1 


COdL 


(90) 


Current  design  technology  takes  this  variation  in  the 
(C n./C.y  natio  into  account  by  "derating"  the  theoretical 
design  value  for  the  number  of  transfer  units  (X©)  depending 
on  the  average  value  of  this  ratio  for  the  matrix.  The 


relationship  is 


29. 


rar** 
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(91) 


X-  X0[i-£c ,/qy‘ "] 


where  \  is  the  number  of  transfer  units 

actually  achieved  in  the  regenerator 
initially  designed  to  provide  a 
value  of  . 

* 

This  equation  is  not  physically  realistic  ,  since 
it  equals  zero  at  a  point  where  ,  and  it  predicts 

negative  values  of  A  if  the  matrix  stream  capacity  is  less  than 
the  gas  stream  capacity.  We  expect  regenerator  performance  to 
suffer  seriously,  of  course,  when  this  ratio  becomes 
small,  but  there  is  no  a  priori  reason  that  heat  exchanger 
performance  should  dramatically  cease  when  the  ratio  is 
unity. 

30 

Lambertson  analyzed  this  variation  of  regenerator 
efficiency  as  a  function  of  stream  capacity  ratio,  assuming 
a  constant  specific  heat  for  the  matrix  lengths,  by 
employing  a  finite  difference  heat  balance  network  and 
iterating  it  using  a  digital  computer.  His  results 
show  conclusively  that  the  number  of  transfer  units  — 
and  hence  the  actually  achieved  efficiency  (£)  --  are 
finite  even  when  the  ratio  is  unity  along  the  entire 
matrix  length.  Note  that  the  actually  achieved  efficiency 
and  the  actually  achieved  value  for  transfer  units  are 
related  by  an  equation  similar  to  the  defining  equation 
for  theoretical  efficiency: 

A 

A- 1  (92) 


4- 


* 

Depending,  of  course,  on  how  the  designer  determines 
the  "average"  values  to  be  assigned  toC^andC^. 
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In  actuality,  since  the  stream  capacity  parameter  of  the 
matrix  varies  along  the  length  of  the  regenerator  due  to 
the  maxtrix  specific  heat  fall-off  with  temperature,  we 
will  find  that  the  ratio  will  (typically)  approach  a  value 
of  unity  only  towards  the  cold  end  of  the  regenerator. 

If  we  consider  the  first  increment  of  gas  flowing 
through  the  regenerator  during  the  "hot  blow"  half  cycle 
(i.e.,  gas  flow  from  toward*"!^  ),  the  gas  will  be  cooled 
so  as  to  approximate  the  roughly  linear  temperature  gradient 
of  the  matrix  until  this  stream  capacity  ratio  attains  some 
low  value  approaching  unity.  At  these  low  values  of  the  ratio 
the  gas  will  continue  to  be  cooled,  but  will  itself  tend  to 
more  rapidly  warm  up  the  lowest  temperature  segment  of  the  matrix 
since  the  specific  heat  fall-off  of  the  solid  now  results 
in  the  gas  stream  capacity  being  comparable  to  —  or  even 


greater  than  —  that  of  the  matrix.  When  the  second  and 
later  increments  of  gas  reach  this  same  point  along  the 
regenerator  length,  the  warmer  matrix  temperature  now  obtain¬ 
ing  means  that  poorer  heat  transfer  will  take  place.  Con¬ 
sequently,  the  overall  effect  results  in  a  regenerator  which 
"looks"  to  be  shorter  than  the  design  length  (since  the 
expected  linear  temperature  gradient  at  the  cold  end  is 
distorted),  so  the  effective  number  of  transfer  units, A  , 
is  less  than  the  theoretical  design  value,  Ae,  as  can  be  seen 
from  Eq.  82  . 

Results  stemming  from  Lambertson's  finite  difference 


network  analysis  can  be  fit  empirically  (for  a  limited 
range  of  A© values  and  )  ratios)  with  an  equation 

of  the  form: 


6-  l-  I  (C%/Ct$  (93) 


Lambertson's  analysis  was  limited  to  conditions  where 
the  (  )  ratio  ranged  from  1-5  and  where  values 


i 

-  W- 
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of  Ko  varied  from  1-10.  Howard  and  Bahnke^"  subsequently- 
used  similar  finite  difference  iteration  techniques  to 
extend  these  results  to  values  of  in  excess  of  200. 

The  range  of  the  Howard  calculations  encompasses  values 
of  the  variables  (  c  *£,  )  and  usually  encountered  in 
cryogenic  refrigerator  design. 

The  form  of  Eq.  93  was  fit  to  the  Howard  data 
set  using  a  multi-linear  regression  technique,  with  the 
resulting  value  for  the  coefficients! 

e*-€„[|-.OS«1  J  (94) 


The  published  data  set  includes  too  few  points  to 
allow  a  comprehensive  assessment  of  the  fit  achieved,  but 
based  on  the  data  available  a  correlation  coefficient  of 
about  95%  was  obtained.  The  empirical  equation  predicts 
efficiency  values  which  differ  by  a  maximum  of  about  ±>5%> 
from  the  Howard  results.  Even  this  small  variance  is 
undesirable,  but  since  experimental  verification  of  the 
published  data  has  apparently  not  been  undertaken  -  at 
least  not  to  this  third  decimal  point  accuracy  -  further 
work  to  improve  the  statistical  fit  is  not  warranted. 

Since  the  value  of  Ct  varies  along  the  length  of  the 
regenerator  it  is  apparent  that  the  actually  achieved 
efficiency  also  varies. 

We  can  express  this  effect  best  with  the  concept  of 
actually  achieved  efficiency  per  segment; 


where 


r  /  u) 

/  /j  i  -) 

J 

(95) 

^  are 

average  values  of  the  efficiency, 

matrix  capacity  and  gas  capacity, 

*  4^ 

respectively,  for  the  A  "  segment 
of  the  regenerator. 


The  smaller  these  (hypothetical)  regenerator  segments, 
the  more  accurate  should  be  the  efficiency  prediction 
derived  from  this  equation. 

For  purposes  of  this  study  program,  we  shall  assume 
that  the  portions  of  the  heat  transfer  supplied  by  the 
first  stage  regenerator  and  the  second  stage  regenerator 
are  idealized,  and  these  two  regenerators  operate  with 
100%  efficiency.  Consequently ,  gas  exiting  the  second 
stage  regenerator  during  the  "hot  blow"  will  be  assumed 
to  be  at  exactly  ,  the  second  stage  operating  temp¬ 
erature ,  as  indicated  in  Figure  30. 

We  shall  arbitrarily  subdivide  the  final  length  of 
this  heat  transfer  gradient,  the  so-called  third  stage 
regenerator,  into  three  segments,  shown  in  Figure  31, 
each  of  which  has  an  efficiency  (  )  associated  with  it; 

each  shall  be  considered  a  function  of  the  capacity 

ratio  only. 

Since  the  mass  flow  through  the  regenerator  is  the 
same  in  both  the  cold  blow  and  the  hot  blow  period,  it 
can  be  shown  that  for  steady  state  operation  the  regener¬ 
ator  efficiency  can  be  alternatively  expressed  in  terms 
of  average  exit  temperatures  instead  of  the  heat  quantities 
of  Eq.  62  .  Thus,  in  general, 


(96) 


or  fpr  the  total  cold  gradient  of  our  model: 

ATJurual 


-=■  (ATV  srfr  *  )  TT*  ST6  +  C&T UTW  * )  j** 

Based  on  thi's  analytical  approach,  the  efficiency  equations 
for  the  three \cold-end  segments  of  Figure  31  can  be  approxi¬ 
mated  as: 


Figure  31.  Idealized  Refrigerator  Model  With  Segmented 
Cold  Gradient  Regenerator 
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(98) 


_  l/*j  —  1/v 


(99) 


*  C£»-T<0 


^3  ~ 


(100) 


— r  Tfrca-T;3 

where  Uv*  is  the  actual  exiting  gas  average  temper¬ 
ature  . 


The  €.3  can  he  computed  from  Eq.  95  sequentially 
for  each  of  the  three  segments  (using  the  segment-average 
for  the  capacity  ratio)  and  the  resulting  value  used  to 
solve  for  the  actual  temperature  difference  produced  in 
each  segment.  This  done,  the  effective  efficiency  for 
the  overall  cold  gradient  regenerator  can  be  computed  using 
the  formula  of  Eq.  97 

where  ~  •  (101) 

Although  this  approach  is  not  rigorously  accurate 
(because  in  actuality  each  of  the  segment  temperatures 
is  time-varying  during  a  "blow”)  it  should  provide  a 
better  prediction  of  regenerator  performance  than  the 
conventional  approaches  of  assigning  a  single  "lumped" 
average  value  to  the  capacity  ratio  based  on  the  expected 
temperature  distribution  of  the  matrix. 
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3*3  Refrigerator  Performance  Calculation 

So  long  as  we  keep  clearly  in  mind  the  physical  model 
just  developed  of  three  refrigerator  stages  operating 
in  parallel,  it  is  convenient  to  treat  the  third-stage 
regenerator  as  if  it  were  a  serial  heat  exchanger,  oper¬ 
ating  between  the  limits  and  lc.  Conventional  design 
techniques  involve  the  determination  of  machine  size 
parameters  based  on  the  efficiency  of  this  third-stage 
regenerator  computed  as: 

which  is  a  logical  extension  of  Eqs.  98,  99  and  100  . 

Heretofore  in  this  report,  all  design  equations 
have  been  derived  based  on  the  refrigeration  produced  as 
the  working  fluid  is  taken  through  one  complete  cycle  of 
the  PV  diagram  of  Figure  26.  in  order  to  apply  these 
equations  to  compute  flow  and  power  requirements  under 
dynamic  operating  conditions,  several  further  approxi¬ 
mations  will  be  required.  We  earlier  defined  the  cycle 
time ,  as  the  time  required  for  the  system  to  achieve 

all  the  states  of  the  P-V  diagram;  consequently,  it  is 
evident  that 

“  -  'he*'*,  -  I  (103) 

where  /Wis  the  number  of  cycles  per  second  and 
N  is  the  RFM  of  the  mechanism. 

It  was  also  noted  that 

do4) 

where ^^is  the  "blow  time"  when  gas  flows 
through  the  regenerator. 


TCa-7, 


ka-Tc 


(102) 
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The  total  mass  of  gas  available  for  expansion  during 
each  cycle  is  then  approximately  the  average  regenerator 
flow  rate  multiplied  by  the  flow  time  or: 


/wv- •AVN'fc^  ■= 

(105) 

The  total  refrigeration  produced  per  second, 
evidently, 

is  then 

^  —  C?  /v^ 

(106) 

Thus ,  from  Eq.  60  , 

(107) 

(108) 

~  ~krrts  &TCS* 

(109) 

'*  “^5.  rr^Cy  • 

(110) 

Regenerator  efficiency  from  Eq. 62  can  likewise  be 
defined  in  terms  of  refrigeration  rate, 

€  -  C  ott  V  >  (ni 

so  that  the  number  of  transfer  units  required  in  order 
for  a  regenerator  to  operate  between  the  limits  Ih  and 
is  found  in  the  following  manner: 
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Since 


(112) 


We  rearrange  and  find, 

.  ^  j _  A»-~^  \ 

W"  JUU&A&  Or  ^  Tc  J 


(113) 


so  long  as  we  recognize  from  Figure  28  that  regenerator 
efficiency  must  exceed  0.99  for  practical  refrigerators. 
The  mass  flow  rate  required  to  provide  the  desired 
refrigeration  is  approximately 


I 

•  Os 

nrr\  —  - - - 

vrvc.  K  CKAO 


(114) 


In  order  to  assess  whether  or  not  small  particle  tech¬ 
nology  can  provide  a  worthwhile  benefit  to  refrigerator 
efficiency,  a  set  of  calculations  has  been  performed 
comparing  the  performance  of  a  conventional  lead  ball 
regenerator  with  a  regenerator  employing  various  sizes  of 
ultrafines.  The  ultrafines  as  pointed  out  in  the  preceding 
section,  achieve  an  improved  heat  capacity'  at  low  temperature 
due  to  the  large  surface  to  volume  ratio  of  the  material. 

If  the  ultrafines  were  merely  packed  into  the  regener¬ 
ator  shell  as  a  powder  bed,  the  pressure  drop  for  required 
values  of  \0  would  be  so  high  that  no  refrigeration  could 

be  produced.  The  equation  to  describe  pressure  drop  through 

22 

a  powder  bed  is  the  Kozeny-Carmen  equation  : 


where 


is  the  volume  specific  surface 

is  the  cross  sectional  area  of  the  bed 

is  the  porosity  of  the  pack 

is  the  length  of  flow  path 

is  the  gas  viscosity 
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rm  is  the  rate  of  flow  corrected  for 

'C 

compressibility 

and  U.  ,  the  Kozeny  constant,  equals  5  for 
spherical  particles. 

It  can  easily  be  verified  that  for  even  small  flow 
values  the  pressure  drop  is  prohibitive.  It  is  evident, 
then,  that  the  ultrafines  must  be  processed  into  bulk 
shapes  to  be  useful  as  regenerator  elements.  The 
feasibility  of  achieving  this  compacting  and  shaping, 
as  well  as  the  effect  this  would  have  on  the  surface  to 
volume  enhanced  specific  heat  contribution  are  discussed 
in  Section  IV. 

For  purposes  of  the  refrigeration  comparison  calcu¬ 
lations,  it  will  be  assumed  that  the  ultrafines  are 
processed  into  a  shape  with  a  hydraulic  radius  comparable 
to  the  conventional  lead  ball  regenerator.  For  0.005  in. 
balls  which  will  be  considered  the  baseline  packing  for 
the  regenerator  matrix  we  find 

51,  -  •  ooi'*<+  c-M.  (116) 


For  comparison,  if  a  fine  mesh  screen  were  used, 

•  OOi 3^  C/v*  (116a) 


For  convenience,  we  will  assume  the  ultrafines  are 
compacted  into  spherical  shapes  of  the  baseline  packing 
size  so  the  enhanced  performance  can  be  directly  compared 
to  a  bulk  lead  regenerator  of  similar  configuration. 

Although  several  materials  were  identified  in 
Section  II  which  evidence  a  greater  surface  specific 
heat  enhancement,  it  is  likely  that  initial  testing  of 
this  effect  will  employ  ultrafines  of  lead  since  its 
toxicity  characteristics  are  well  known  and  there  is 
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some  experience  to  be  drawn  upon  in  both  handling  the 
material  and  producing  powders  of  the  required  size. 

Thus,  the  curve  of  Figure  25  can  be  used  to  determine 

<? 

the  specific  heat  characteristics  of  the  heat  exchanger. 

These  values  are  tabulated  in  Table  3  for  the  four 

test  cases  for  which  refrigerator  performance  calculations 

have  been  performed: 

Case  I  -  T„  =  35°K,  Tc  =  15°K 

Case  II  -  TH  =  35°K»  Tq  *  10°K 

Case  III  -  TH  «=  35°K,  Tc  *  5°K 

Case  IV  -  T„  **  15°K,  T„  =  5°K, 

It  c 

In  order  to  obtain  numerical  results  for  refrigerator 
performance,  the  baseline  model  of  Table  4  was  assumed. 

The  values  in  the  table  are  intended  to  be  compatible 
with  calculations  for  the  High  Capacity  Vuilleumier 
Refrigerator  performance,  as  determined  by  the  Vuilleumier 
Cycle  Cryogenic  Refrigerator  Optimization  Computer 
Program^2  That  program  is  tailored  to  the  Hi-Cap  VM 

machine  whereas  the  calculations  reported  here  are  based 
on  the  idealized  Stirling  cycle  of  Figure  26,  so  power 
estimates  and  sizing  can't  be  directly  compared  between 
the  two.  However,  the  Hi  Gap  results  allow  us  to  assign 
physically  realistic  values  to  our  model  for  gross  heat 
loads  and  operating  temperatures  for  the  three  stages, 

In  particular,  it  should  be  noted  that  the  ratio  of 
regenerator  heat  loss  to  gross  third  stage  refrigeration, 
f  determined  by  that  program  and  as  used  in 

Table  4  is  typical  of  well  designed  Stirling  eycle  or 
VM  cycle  machines,  and  the  gross  third-stage  refrigeration 
of  2.785w  is  a  typical  value  needed  to  produce  a  net  of 
about  0.3w  available  for  the  external  load  at  about  10°K. 
The  difference  between  the  gross  output  and  net  output  is 
the  sum  of  the  regenerator  loss,  whieh  we  primarily  focus 
on  in  this  report,  and  the  other  stage  losses  treated  in 
detail  elsewhere^. 
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Table  3.  Specific  Heat  Values  for  Third  Stage  Regenerator 

(Tabulated  Values,  in  Joules/fc.Ma/K,  Correspond  to 
Endpoints  of  Segmented  Model  —  Fig.  23) 


TABLE  4 

BASELINE  REFRIGERATOR  MODEL 


1)  Crankcase  temperature 

2)  First-stage  temperature 

3)  First-stage  gross  output 

4)  Second-stage  temperature 

5)  Second-stage  gross  output 

6)  Regenerator  efficiency  for  first  stage 
and  second  stage 

7)  Mechanism  operating  speed 

8)  Base  cycle  pressure 

9)  Expansion  pressure  ratio 

10)  Third-stage  parameters. 

a)  Regenerator  hydraulic  radius  (.0 05  lead 

spheres) 

h)  Regenerator  Reynolds  Number  (optimized  so 

f/j  is  at  16.6) 

0 

c)  Q3-.4  (gross  output) 

e)  Net  refrigerator  output 

f)  Maximum  regenerator  pressure  drop 


324°K 

67  °K 
57  w 
35°K, 

Case  1,  2,  3 
15°K  Case  4 

22  w 

100 % 

600  RFM 
40  atm 
1.4 

0.0013  cm 
500 

2.785  w 
0.395 
0.3  w 
0.5  atm 
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3.3*1  Interpretation  of  Calculated  Results 

The  results  of  refrigerator  performance  analysis 
are  presented  in  Table  5*  For  each  test  case,  the 
required  third-stage  values  of /Yr^^and  to  satisfy 

the  baseline  model  conditions  of  Table  4  were  computed 
from  Eqs,  113  and  114  . 

Using  an  iterative  approach  a  design  value  for  heat 
transfer  units  (  \0)  for  the  (bulk  material)  lead  sphere 
matrix  was  assumed,  then  sizing  calculations  were  performed 
and  the  appropriate  derating  factors  (to  account  for  low 
temperature  specific  heat  fall-off)  were  applied.  The 
resulting  value  of  "effective"  number  of  transfer  units 
was  compared  with  the  value  (from  Eq.  114  ,  and 

the  iteration  was  continued  until  a  reasonable  agreement 
was  achieved.  This  iterative  technique  determined  the 
"design"  value,  ,  to  be  employed  for  sizing  the  physical 
parameters  of  the  mechanism. 

For  each  test  case  three  additional  iterative 
computations  were  performed!  one  each  for  the  spheres 
made  of  100A°,  50A°  and  22A°  lead  ultrafines.  The  iterative 
technique  again  showed  an  optimum  design  value  associated 
with  that  size  of  ultrafines,  to  be  employed  for  sizing  the 
mechanism.  The  optimized  values  of  give  rise  to  values 
for  effective  efficiency  which  differ  between  the  bulk 
material  spheres  and  the  compacted  spheres  of  ultrafines. 

This  incremental  change  in  efficiency  is  the  basis  for  a 
potential  savings  in  input  power,  for  a  reduction  in  physical  si 
or  for  an  increase  in  the  refrigeration  available  to  the 
external  heat  load  --  thus  providing  a  kind  of  safety 
faotor  to  allow  for  better  temperature  control  of  the  load. 

It  is  important  to  understand  tie  implications  of  the 
increased  efficiency  and  the  options  for  improvement  which 
are  presented, 
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Option  1.  Reduction  in  the  Design  Value  of  Transfer  Units 
In  general,  an  enhanced  specific  heat  obtained  through 
the  use  of  ultrafines  in  the  regenerator  matrix  allows 
the  refrigerator  designer  to  select  a  lower  design 
value  for  than  if  he  were  working  with  conventional 
matrix  materials.  Several  parameters  can  be  traded 
off  to  optimize  the  method  of  lowering  the  design  value 
but  Eq.  117  gives  some  insight  as  to  the  benefit 
which  will  accrue. 


I  »  Wy* 

J 


(117) 


The  variables  are  to  some  extent  functionally  inter¬ 
related,  but  it  can  be  seen  that  a  lower  value  of 
allows  the  designer  to  either  proportionally  reduce 
the  regenerator  length  or  to  increase  the  hydraulic 
radius.  In  the  one  case  he  is  opting  to  reduce  the 
size  of  his  machine  and  in  the  other  to  reduce  the 
flow  friction  losses  in  the  regenerator.  In  addition 
to  reducing  the  thermodynamic  flow  friction  losses 
this  latter  approach  will  contribute  to  improved 
reliability  by  lowering  the  loads  on  bearings  and 
seals.  This  first  option  is  not  considered  to  be 
the  most  advantageous  way  of  capitalizing  on  the 
specific  heat  enhancement. 


Option  2.  More  Capacity  at  T  for  the  External  Load 

For  a  given  design  value  for  Xo  »  an  improved  matrix 
specific  heat  results  in  a  greater  value  for  regen¬ 
erator  efficiency »  since  the  effective  value  for  X 
increases.  If  the  regenerator  efficiency  is  improved, 
then  the  regenerator  losses  are  reduced, 

Ob~-  0"O  Qvi  (US) 
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This  reduction  in  regenerator  heat  loss  is  directly 
available  as  added  capacity  for  the  external  load 
since  from  Eq.  66  : 

Qxt*  (119) 

For  the  same  gross  refrigeration,  and  no  change  in 
stage  losses,  a  reduction  in  Q*  directly  increases 
the  available  Qo+-  If  the  designer  elects  this 
option,  he  will  still  design  to  the  same  X.  value, 
but  will  achieve  a  safety  factor  over  the  designed 
net  capacity  as  inuicated  in  the  table. 


Option  3*  Reduction  of  Input  Power  for  the  Same  Safety  Factor 
in  Net  Refrigeration 

If  the  designer  of  current  technology  refrigerators  must 
provide  the  same  safety  factor  which  could  be  achieved 
automatically  through  the  use  of  ultrafine  particles  in 
the  regenerator  matrix,  he  must  increase  the  gross 
refrigeration  to  overcome  the  poorer  efficiency  of 
conventional  materials.  The  increase  in  -gross  output 
must  be  achieved  through  increased  mass  flow  rate,  which 
in  turn  produces  a  larger  absolute  value  of  regenerator 
loss  for  the  same  effective  efficiency  value.  Eq.  70 
can  be  rearranged  as, 


VR= 


Op 


*****  \ 


(120) 


so  that  the  power  requirement  can  be  computed  from  the 
regenerator  loss  value  and  the  effective  number  of  transfer 
units.  Values  of  input  power  are  compared  in  the  table 
versus  a  percentage  safety  factor  provided,  a)  automatically 
through  the  use  of  ultrafines,  and  b)  by  way  of  increased 
mass  flow  when  conventional  technology  is  employed. 
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Table  5.  Predicted  Regenerator  Improvements  Obtainable 
From  Small  Particle  Technology,  (sheet  1  of  2) 


Using  Test  Case  II  of  Table  5  to  illustrate  these 
options,  the  following  conclusions  can  be  drawn  from 
values  tabulated. 

o 

1.  If  22  A  particles  are  used  to  form  the  .005  in 

diameter  lead  spheres  in  the  regenerator  matrix, 

the  design  value  of  A©  for  this  35°K  to  10°  stage 

can  be  reduced  from  167  (required  for  solid  lead 

spheres)  down  to  131«  This  reduction  shortens  the 

regenerator  from  4.34  cm  down  to  3*40  cm,  decreases 

2  2 

the  critical  flow  area  from  .117  cm  to  .103  cm  , 
and  reduces  the  lead  volume  requirements  in  the 
matrix  from  .83  cm^  to  .57  cm^. 

2.  Alternatively,  if  the  original  sizing  calculations 
are  used  and  the  refrigerator  uses  a  design  value 
for  A0  of  I67,  the  spheres  made  of  22  £  particles 
will  provide  an  effective  value  for  transfer  units 
of  66.7,  as  compared  with  the  bulk  material  value 
of  38.  Additionally,  the  regenerator  loss  will 
decrease  to  .6l80Vu.,  as  opposed  to  1.0713  W-  for  the 
bulk  material,  providing  a  151$  safety  factor  for 
the  0.3w.  net  refrigeration  requirement. 

3.  Using  this  same  design  value  of  167  for  A®  ,  the 
151$  safety  factor  will  be  achieved  with  a  refrig¬ 
erator  power  input  of  471w*  •  If  the  bulk  material 
were  used,  and  a  comparable  safety  factor  were 
desired,  mass  flow  rate  would  have  to  be  increased 
and  the  minimum  input  power  necessary  to  achieve 
the  needed  gross  refrigeration  increase  would  be 
494\m  . 
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IV 


ULTRAFINE  PARTICLE  TECHNOLOGY 


4.1  Introduction 

In  an  earlier  section  of  this  report  it  was  deter¬ 
mined  that  finely  divided  powders  show  an  increase  in  the 
low  temperature  specific  heat  value  when  compared  with 
the  bulk  material.  This  effect  is  at  least  approximately 
proportional  to  specific  surface  of  the  sample.  This 

proportionality  probably  holds  at  least  down  to  particle 

o 

diameters  of  20  -  30  A. 

For  maximum  improvement  in  the  performance  of  cryogenic 
refrigerators  it  is  obviously  desirable  to  maximize  the 
specific  heat  of  the  material  used  to  store  energy  in 
the  regenerator  matrix.  Thus  it  is  of  interest  to  deter¬ 
mine  whether  this  desired  improvement  can  be  effected  by 
utilizing  small  particle  technology  to  take  advantage  of 
this  specific  heat  enhancement. 

Calculations  presented  earlier  indicate  that  a  worth¬ 
while  improvement  in  regenerator  performance  due  to  this 

0 

effect  would  require  that  particles  100  A  in  diameter  or 

smaller  be  used  to  make  up  the  heat  exchanger  matrix. 

Particles  of  this  size  category  are  usually  referred  to  as 

ultrafines,  and  are  comparable  in  size,  at  least  in  the 

low  end  of  the  range,  to  some  of  the  larger  single  molecules. 

Consequently,  we  are  approaching  a  physical  limit  of  what 

can  be  still  called  a  particle  when  we  consider  substances 
o 

much  below  100  A. 

In  the  past  decade  numerous  technical  reports  have 
been  published  regarding  the  production,  measurement,  and 
handling  of  small  particles  of  both  dielectric  and  metal 
substances.  Because  metal  powders  have  potentially  a 
larger  specific  heat  for  a  given  Debye  temperature  value, 
generally  higher  bulk  densities,  and  generally  more  favorable 
mechanical  properties,  they  are  more  attractive  than  dielectrics 
for  use  in  low  temperature  regenerator  experiments. 

Much  of  the  early  work  in  production  and  handling 
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techniques  for  ultrafine  powder,  however,  involved  dielectric 
materials  instead  of  metals.  These  materials  are  not 
generally  as  hazardous  to  work  with  as  the  pure  metals, 
which  tend  to  be  highly  reactive  as  well  as  being,  in  many 
cases,  highly  toxic. 

4.2  Production  of  Particles 

There  are  a  surprisingly  large  number  of  ways  that 

ultrafine  particles  can  be  produced.  Rather  than  survey 

all  of  the  possible  methods  in  this  report,  the  reader  is 

cited  to  several  excellent  review  articles  which  themselves 

35 

provide  numerous  references 

Several  of  the  methods,  including  the  electrostatic 
precipitation  and  inert  gas  evaporation  methods,  seem  to 
show  up  repeatedly  in  the  literature.  We  can  thus  con¬ 
clude  that  the  resulting  particle  production,  especially 
from  the  inert  gas  technique  is  more  or  less  reliable  and 
that  the  size  distribution  is  narrow  enough  to  satisfy 
scientific  measurement  requirements. 

Most  of  the  ultrafine  powder  samples  first  used  for 
specific  heat  measurements  of  dielectric  materials  in  the 
1950s  were  produced  by  an  electrostatic  precipitation 
process,  wherein  the  crystal  was  heated  in  a  dry  atmos¬ 
phere  until  it  vaporized.  The  smoke  particles  thus 
produced  were  collected  in  a  high  voltage  alternating 
current  electrostatic  precipitator.  They  were  subsequently 
handled  under  scrupulously  anhydrous  conditions  since 
exposure  to  moisture  caused  an  increase  in  particle  size. 

Specific  surfaces  of  sodium  chloride  as  large  as  96  m  /gm 

0 

(corresponding  to  about  1000  A  cube  edge)  have  been  pre- 
pared  by  this  method.  Smaller  particles  of  MgO,  having  a 
specific  surface  of  166  m  /gm  (corresponding  to  about  100 
cube  edge),  were  produced  by  Lien  and  Phillips'*  using  a 
dehydration  technique,  the  details  of  which  were  not 
reported  by  them. 

Other  methods  which  reportedly  yield  predictable 
and  reproducible  results  include  the  decomposition  of 
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specially  prepared  metal-organic  compounds^?  —  apparently 
used  primarily  for  the  preparation  of  refractory  oxides  -- 
and  the  high  temperature  decomposition  of  carbonyl  vapor^® , 
which  was  an  early  method  used  to  produce  particulates  of 
high  purity  metal. 

Much  of  the  recent  work  with  metals  involves  one 
variation  or  another  of  the  inert  gas  evaporation  method. 

In  this  technique  a  pure  metal  sample  is  homogeneously 
and  rapidly  heated,  often  either  by  induction  heating  or 
temperature-stabilized  oven  techniques,  to  produce  a  super¬ 
saturated  vapor  within  an  inert  gas  atmosphere  of  the 
evaporation  apparatus. 

Nucleation  and  growth  of  the  metal  particles  begins 
in  the  gas  atmosphere.  During  this  process  traces  of 
impurities  are  highly  undesirable  since  they  will  cause 
oxides,  nitrides,  or  other  contaminated  forms  of  the 
particles  to  grow.  Success  of  this  method  in  producing 
narrow  size  distributions  of  ultrafine  high  purity  particles 
seems  to  be  related  to  the  size  of  the  inert  gas  molecule, 
as  well  as  to  the  ability  to  rapidly  achieve  a  highly  super¬ 
saturated  condition  of  the  vapor.  This  high  degree  of 
supersaturation  is  conducive  to  the  formation  of  a  large 
number  of  condensation  nucleii  and  to  subsequent  slow 
growth,  which  allows  a  reasonable  (but  short)  time  to 
collect  the  ultrafines  --  often  on  a  water-cooled  surface 
which  inhibits  further  nucleii  development.  The  particles 
so  collected  must  be  stored  in  vacuum  or  in  an  inert  gas 
atmosphere  since  the  high  surface  to  volume  ratio  results 
in  a  large  excess  surface  energy  which  not  only  tends  to 
form  agglomerations  of  particles  but  also  results  generally 
in  strong  pyrophoric  tendencies. 

An  interesting  variation  on  this  inert  gas  evaporation 
technique  will  be  discussed  in  a  later  section  as  the  best 
candidate  system  for  providing  a  production  scale-up  of 
the  technology,  should  the  later  test  phases  of  this  program 
confirm  the  current  analytical  predictions  of  refrigerator 
efficiency  improvement. 


4.2.1  Measurement  and  Grading  of  Ultrafine  Powders.  The  BET 
adsorption  isotherm  technique,  introduced  earlier,  has 
been  the  accepted  standard  method  of  determining  the  area 
of  finely  divided  powders  since  its  inception.  In  this 
method  the  amount  of  gas  (typically  nitrogen  maintained 
near  its  condensation  point)  which  adsorbs  on  the  active 
surface  of  the  powder  is  determined  by  measuring  a  change 
in  pressure  of  the  test  chamber  before  and  after  the 
adsorption  process  takes  place.  This  method  is  time  con¬ 
suming  and  requires  precise  process  controls,  but  provides 
a  highly  accurate  result. 

There  have  been  recent  attempts  to  develop  less  time 
consuming  dynamic  gas  flow  techniques  utilizing  air  per¬ 
meability  measurements  to  determine  powder  surface  areas, 
but  in  general  the  accuracy  is  questionable  and  the  range 
of  measurement  capability  doesn't  extend  to  the  ultrafines 
which  are  of  interest  for  this  program. 

An  important  disadvantage  of  the  BET  method  is  that  it 

measures  only  total  surface  area,  and  doesn’t  allow  a 

determination  of  the  distribution  of  particle  sizes  within 

a  sample.  Size  distribution  is  critically  important  for 

the  regenerator  specific  heat  enhancement  application 

.  o 

because  only  those  particles  below  about  100  A  in  size 

will  evidence  a  large  enough  enhancement  effect  so  as  to 

improve  regenerator  efficiency.  Consequently,  if  a  given 

test  sample  of  suitable  "average”  particle  size  were  made 

up  of  an  extreme  bi-modal  distribution,  e.g.  one  group  of 
o  o 

300  A  and  one  of  10  A  particles,  the  resulting  measurable 
enhancement  in  specific  heat  could  conceivably  be  almost 
nil ,  depending  on  the  numbers  of  members  of  each  distri¬ 
bution  which  combined  to  produce  the  average  surface  area 
of  the  sample. 

Because  of  the  growing  industrial  interest  in  ultra- 
fine  powders  for  high  performance  magnets,  self-lubricated 
bearings,  and  dispersion  strengthened  alloys,  considerable 
experimental  and  theoretical  work  in  the  past  decade  has 
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addressed  the  question  of  size  distribution. 

Since  the  BET  method  cannot  investigate  the  size 
distribution  function,  electron  microscope  techniques  are 
commonly  employed  for  this  purpose.  Although  many  re¬ 
searchers  also  employ  this  instrument  to  estimate  the  actual 
specific  surface  area  of  a  powder  sample,  chere  is  apparently 
no  consensus  that  it  can  replace  the  BET  method  for  area 
determination. 

Granqvist  and  Buhrman^  have  reported  that  ultrafine 

powders  tend  to  take  a  spherical  shape,  rather  than 

follow  the  appropriate  crystal  habit,  if  the  particle 

o 

diameter  is  below  about  200  A.  They  further  determined  that, 
provided  the  inert  gas  evaporation  technique  is  used, 
the  size  distribution  follows  a  log  normal  distribution 
function  with  standard  deviation  in  the  range  1.36^<T^  1.60. 

They  corroborated  this  theory  with  extensive  electron 
microscope  studies  of  ultrafines  produced  from  aluminium 
oxides,  chromium,  iron,  cobalt,  nickel,  copper,  gallium, 
magnesium,  zinc,  and  tin. 

Consequently,  if  the  inert  gas  evaporation  method  of 
producing  ultrafines  is  employed,  the  question  of  how  to 
grade  and  separate  the  particles  to  be  used  in  later  test 
phases  of  the  current  program  disappears. 

4.2.2  Predicted  Availability  of  Test  Materials.  There  is  no 
longer  any  question  that  gram-size  samples  of  spherical 
ultrafine  metal  particles,  with  nominal  diameters  of  about 
50  X,  and  acceptably  narrow  size  distributions,  can  be 
produced.  In  order  to  advance  into  the  rov+  phase  of  this 
program,  however,  it  is  necessary  to  determine  how  to  acquire 
such  powders  made  from  the  aterials  identified  earlier  as 
the  best  candidates  for  the  refrigerator  efficiency  improve¬ 
ment  application. 

The  powder  metal  industry  does  not  currently  work 
with  powders  much  smaller  than  400  mesh.  A  400  mesh  screen 
has  openings  of  about  4 5Unw  on  a  side.  Particles  which 
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pass  through  such  a  screen  are  obviously  smaller  than 
the  opening,  but  still  two-to-three  orders  of  magnitude 
larger  on  average  than  the  ultrafines  which  will  exhibit 
a  sufficiently  enhanced  specific  heat  effect  for  purposes 
of  the  regenerator  application. 

An  extensive  survey  of  the  leaders  in  this  powder  metal 
industry  has  demonstrated  convincingly  that  50  -  100  X 
powders  are  not  currently  available  off  the  shelf, 

A  survey  of  industrial  research  laboratories  which 
specialize  in  materials  development  also  turned  up  negative 
results.  This  survey  included  companies  like  Union  Carbide, 
Celanese,  Dupont,  Bethlehem  Steel,  Hercules  Inc.,  GCA  Corp. , 
and  Ethyl  Corp. ,  among  others.  Of  this  group  only  Ethyl 
Corp.  reported  prior  experience  in  ultrafine  particle 
production  (of  aluminum).  Their  activity  was  discontinued 
some  time  back  for  economic  reasons. 

Finally,  a  survey  of  academic  institutions  was  con¬ 
ducted,  using  author  names  and  institutions  obtained  from 
recently  published  papers  on  small  particles. 

Based  on  this  effort  it  was  determined  that  Cornell 
University,  Ithaca,  N. Y. ,  has  the  capability  currently  to 
produce  particles  of  the  size  and  distribution  character¬ 
istics  necessary  for  further  experimental  work  on  this 
program.  Based  on  discussions  with  personnel  from  that 
organization,  indium  must  be  eliminated  from  the  list  of 
candidate  materials  presented  earlier  since  its  ultra- 
fines  are  apparently  exceedingly  soft  and  tend  to  form 
films  instead  of  collectible  particles.  Lead  and  selenium, 
our  other  two  primary  candidates,  are  considered  by  them 
to  be  viable  materials  for  particle  production,  with  the 
toxicity  characteristics  of  the  latter  needing  more 
investigation. 

As  noted  earlier,  the  preferred  method  of  particle 
preparation  is  the  inert  gas  evaporation  technique  which 
was  analyzed  in  detail  by  the  researchers  at  Cornell. 

For  production  scale-up  of  the  particle  production, 


109 


providing  that  test  results  in  later  program  phases 
remain  encouraging,  a  variation  of  the  inert  gas  technique 

has  already  been  pioneered  by  a  Japanese  laboratory. 

40 

Yatsuya,  Mihama,  and  Uyeda  recently  reported  on  a 
vacuum  evaporation  technique  which  captures  the  newly 
formed  ultrafines  in  the  oil  film  flowing  outboard  under 
the  influence  of  centrifugal  force  on  a  rotating  disc 
inside  the  bell  jar.  This  is  essentially  a  continuous 
production  operation,  and  the  particles  are  subsequently 
separated  from  the  oil  by  distillation  techniques.  The 
researchers  have  projected  that  the  technique  is  suitable 
for  reasonably  large  scale  production. 

Methods  of  Regenerator  Production 

Very  little  work  has  been  done  to  date  to  develop 
techniques  for  handling  ultrafine  powders  and  forming 
or  shaping  materials  made  from  them.  Since,  as  noted, 
the  particles  can't  be  employed  in  a  refrigerator  merely 
as  a  packed  power  bed  due  to  pressure  drop  considerations, 
it  is  necessary  to  fabricate  special  regenerator  elements 
from  them. 

In  the  process  of  forming  the  particles  into  usable 
shapes,  care  must  be  taken  not  to  destroy  the  high 
surface  to  volume  ratio  which  gives  rise  to  the  enhanced 
specific  heat  effect. 

Several  classes  of  methods  suggest  themselves  for 
forming  the  ultrafines  into  usable  configurations.  First, 
techniques  from  the  powder  metallurgy  industry  can  be 
applied.  In  forming  powdered  metal  components,  a  two- 
stage  approach  is  commonly  followed.  The  raw  powders 
are  initially  pressed  into  a  "green  compact”  stage, 
where  the  amount  of  surface  welding  or  alloying  is  minimal 
and  is  related  to  the  pressure  of  compacting.  This  green 
compact  has  generally  poor  mechanical  properties  until 
it  is  heated  to  just  under  the  melting  temperature  and 
"sintered",  wherein  heavier  surface  welding  takes  place. 


Heavy  surface  welding  is  obviously  undesirable 
for  our  application,  since  discrete  particle  character¬ 
istics  must  be  preserved.  Thus,  an  optimum  green  compacting 
pressure  must  be  determined  to  provide  maximum  strength 
of  the  element  with  an  acceptable  level  of  alloying  or 
welding  as  determined  from  specific  heat  measurements. 

There  is  a  precedent  for  this  approach.  The 
experiment  of  Comsa,  et  al,  on  palladium  ultrafines, 
noted  earlier,  was  designed  so  that  the  powders  were 
compacted  into  a  small  pellet  before  the  specific  heat 
was  measured.  This  processing  yielded  a  Montroll  effect 
which  fit  the  analytical  predictions  exceedingly  well. 
However,  they  did  not  make  any  determinations  of  the 
mechanical  strength  of  the  compact. 

An  advantage  in  this  approach  is  that  the  metal 
compact  will  act  as  a  helium  sponge  in  the  regenerator, 
since  the  large  surface  to  volume  ratio  will  encourage 
an  adsorbent-adsorbate  relationship,  A  disadvantage  in 
the  approach  is  that  there  will  be  a  tendency,  the 
seriousness  presently  unknown,  for  the  lightly  welded 
ultrafines  to  break  off  from  the  compact  and  move  with 
the  gas  flow  through  the  regenerator  as  a  smoke,  in 
addition  to  causing  an  undersirable  virtual  heat  leak 
through  this  "pumping"  action,  the  ultrafines  conceivably 
could  cause  bearing  and  seal  deterioration  if  they  are 
not  contained  inside  the  regenerator. 

With  present  technology,  it  is  impossible  to 
positively  contain  free  particles  smaller  than  about 
200  1.  Consequently  if  the  individual  particles  break 
away  from  the  compact  some  migration  will  take  place, 

One  interesting  approach  to  this  problem  comes  from 
the  field  of  powder  metallurgy  itself,  itudieg  into  the 
porosity  of  metal  oompaots  employ  a  process  called 
mercury  porosimetry,  In  this  technique,  liquid  mercury 
under  high  pressure  is  forced  into  the  pore  structure 
of  the  formed  element,  The  pore  volume  is  determined 
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by  the  amount  of  mercury  which  flows  during  the  pressure 
cycle.  Since  the  pore  structure  is  the  mirror  image  of 
the  compact,  we  can  expect  that  the  mercury  will  flow 
into  a  three  dimensional  shape  which  itself  is  character¬ 
ized  by  a  high  surface  to  volume  ratio.  Current  technology 
allows  mercury  to  penetrate  pore  sizes  as  small  as  18  X 
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with  pressures  achievable  with  commercial  equipment 

Prom  Figure  17  it  is  known  that  mercury  has  a  better 
low  temperature  specific  heat  than  lead,  and  also  that 
the  surface  enhancement  for  large  surface  to  volume  ratios 
is  the  highest  of  any  material  surveyed  to  date.  Con¬ 
sequently,  if  the  porosimetry  approach  does  indeed  yield 
the  mirror  image  structure,  the  potential  for  a  compli¬ 
mentary  effect  to  the  lead  enhancement  is  available. 

Since  mercury  easily  alloys  with  many  metals,  includ¬ 
ing  lead,  it  can  be  expected  that  surface  alloying  will 
take  place  when  the  two  metals  are  in  contact.  This  will 
increase  the  strength  of  the  green  compact,  and  probably 
eliminate  any  objections  to  using  mercury  in  the  regenerator 
applications  founded  on  its  vapor  pressure  characteristics. 

What  is  currently  an  unknown,  however,  is  the  effect 
this  surface  alloying  will  have  on  the  specific  heat 
enhancement  effect.  It  is  believed  that  the  alloy  inter¬ 
face  will  form  a  wave  reflection  or  scrambling  boundary 
for  the  low  frequency  acoustic  vibrations,  and  that  this 
will  allow  the  full  Montroll  effect  to  be  observed.  Until 
the  tests  are  conducted,  however,  this  is  speculation. 

If  the  theory  developed  in  Section  II  for  the  cluster 
structure  materials  proves  accurate,  however,  those 
materials  provide  some  evidence  that  the  surface  enhance¬ 
ment  effect  survives  when  ultrafine  particles  are  embedded, 
and  even  alloyed,  into  a  metal  matrix  binder. 

The  next  most  feasible  approach  to  containing  the 
particles  inside  a  regenerator  draws  on  the  current 
state  of  thin  film  membrane  technology.  A  survey  of 
several  organizations  involved  in  both  cellulose  acetate 


and  teflon  membrane  research  has  provided  encouragement 
that  a  quantity  of  ultrafines  can  be  mixed  into  a  batch 
of  the  membrane  substance  and  rolled  into  thin  sheets. 
Estimates  of  achievable  particle -to-membrance  density 
in  the  50  ~  80%  range  (by  volume)  have  been  obtained. 

There  are  several  possible  ways  of  utilizing  the 
sheet  material.  Membranes  of  this  type  are  commonly 
used  as  filters  or  separators.  They  can  perform  this 
function  because  a  network  of  precisely  sized  pores 
forms  when  the  substance  cures  initially.  The  pore 
size  can  be  controlled  over  wide  limits,  with  the 
smallest  diameter  being  about  0.2^UY»*for  the  cellulose 
acetate  and  .035^*  f  or  the  teflon.  Since  400  mesh  screen 
has  a  pore  size  in  the  tens  of  microns  and  a  hydraulic 
radius  comparable  to  the  lead  spheres  currently  used  in 
regenerators,  a  membrane  formed  to  approximate  a  screen 
of  this  size  could  potentially  serve  as  a  regenerator 
element  "wafer" . 

Advantages  to  this  permeable  membrane  approach 
include  ease  of  forming  the  regenerator  elements  and 
positive  containment  of  the  particles.  Disadvantages 
include  the  fact  that  mechanical  strength  of  a  membrane 
is  low,  and  pressure  cycling  of  the  helium  gas  in  a 
regenerator  could  pose  a  problem  over  time. 

Lastly,  it  is  attractive  to  consider  merely  "gluing" 
the  ultrafines  together  and  avoiding  the  problems  and 
unknowns  inherent  in  the  powder  metallurgy  and  thin 
membrane  technology  approaches  just  discussed. 

It  has  been  suggested  that  an  organic  binder  like 
GE7031  varnish,  which  has  a  good  thermal  conductivity  at 
low  temperature,  could  be  treated  chemically  so  that  it 
would  form  into  a  hard  compact  upon  curing.  Using 
centrifuge  techniques,  the  density  of  ultrafines  could 
be  raised  to  a  sufficiently  high  level  in  the  liquid 
varnish,  and  then  droplets  of  approximately  the  desired 
size  for  regenerator  elements  could  be  formed  and  cured. 


In  addition  to  the  GE  varnish  approach,  a  low  melting 
point  ceramic  which  itself  has  favorable  low  temperature 
specific  heat  characteristics  has  been  suggested  as  such 
a  binder  . 

Of  the  various  methods  discussed,  only  the  green 
compact  approach  has  any  experimental  evidence  to  directly 
support  its  viability.  Consequently,  based  on  work  by 
Comsa,  et  al  noted  earlier,  the  best  candidate  for  form¬ 
ing  the  ultrafines  into  suitable  regenerator  elements 
is  this  one. 

It  is  unavoidable  that  a  significant  amount  of 
developmental  work  must  be  invested  in  this  aspect  of 
the  application  for  the  program  to  progress  beyond  this 
first  phase. 

4.4  Safety  Considerations 

The  pyrophoric  tendencies  of  ultrafine  metal  powders 
have  been  noted  earlier.  In  addition  to  this  potential 
fire  and  explosion  danger,  however,  the  inhalation  of 
these  powders  can  produce  a  significant  health  hazard. 

For  obvious  reasons,  most  inhalation  studies  have 
focussed  on  particles  larger  than  the  ultrafines  of 
interest  fop  regenerator  applications,  Although  mono- 
disperse  particles  in  the  20  -  10Q  &  size  can  now  be 
reliably  produced,  fractional  micron  particles  are 
more  commonly  available  and  have  been  used  extensively 
for  clinical  experiments. 

Based  on  the  work  done  to  date,  the  number  of 
inhaled  partioles  deposited  in  the  respiratory  tract 
reaches  a  minimum  value  at  equivalent  sphere  diameters 
of  ,  The  deposition  fraction  increases  on 

both  sides  of  this  value,  implying  that  in  the  ultra- 
fine  size  region  the  inhalation  problem  must  be  regarded 
as  a  danger, 

Since  lead  is  likely  to  be  the  initial  material 
tested  in  later  phases  of  this  program  some  discussion 
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of  its  hazards  is  warranted.  Lead  can  enter  the  body 
via  the  skin,  by  ingestion,  or  by  inhalation.  Since 
inhalation  is  the  most  common  cause  of  lead  poisoning, 
and  lead  vapors  are  probably  comparable  in  particle 
size  to  the  ultrafines  to  be  employed  on  this  program, 
the  Occupational  Safety  and  Health  Administration  guide¬ 
lines  for  atmospheric  content  of  lead  can  serve  as  a 
starting  point  for  program  safety  considerations.  Any 
exposure  of  humans  to  concentrations  of  lead  exceeding 
0.6  milligrams  per  cubic  meter  of  air  is  considered  a 
serious  health  hazard,  and  the  highest  permissible 
average  concentration  over  an  eight  hour  day  is  0.2 
milligram  per  cubic  meter* 

Safety  procedures  must  also  be  implemented  to 
periodically  monitor  the  lead  content  of  blood  and  urine 
samples  from  all  personnel  working  in  the  vicinity. 

Lead  is  not  on  the  State  of  California,  Division  of 
Industrial  Safety,  list  of  carcinogenic  substances,  so 
no  special  licensing  requirements  must  be  met  to  work 
with  the  substance.  Prudent  experiment  design,  however, 
will  employ  vented  glove  boxes  and  approved  personal 
respirators  for  personnel  in  direct  contact  with  the 
experiment  itself. 


1 


ih’dCLIiSlONS 

Small  Pa rticle  Specific  Heat 

Although  early  measurement  programs  on  the  low 

tempera  lure  specific  heat  of  finely  divided  materials 

P'ovj.dea  i  neons  is  tent  results  and  were  plagued  with 

experimental  difficulties ,  recent  work  in  the  United 

Si.at.es  and  abroad  confirms  that  a  significant  specific 

tea  enhancement  is  available  for  ultrafines  with  an 

o 

equivalent  diameter  below  100  A. 

theoretical  work  to  predict  this  enhancement  has 
neon  under  development  for  nearly  three  decades. 

Rigorous  models  from  lattice  dynamics  technology  are 
complex  to  use  and  currently  do  r.o  better  at  predicting 
the  magnitude  of  the  specific  heat  effect  than  simple 
models  patterned  after  conventional  Debye  analysis. 

Based  on  its  excellent  correlation  to  experimental 
specific  heat  values  for  palladium  ultrafines,  the 
modified  Montroll  equation  presented  in  Section  II  is  an 
accurate  enough  predictor  of  enhanced  specific  heat 
values  i'or  purposes  of  this  intended  application. 

Applying  that  equation  to  the  group  of  low  Debye 
temperature  substances  which  are  candidates  for  use  in 
later  experimental  phases  of  this  program,  we  conclude 
that  lead  is  the  best  choice  for  initial  experimental 
work . 

email  Particle  Availability 

There  is  no  doubt  that  ultrafine  particles  of 
o 

nominally  50  A  diameter,  with  an  acceptably  narrow 
size  distribution,  can  be  produced  using  the  inert  gas 
evaporation  technique.  Since  the  powdered  metal  industry 
is  not  yet  working  with  sub-micron  particles  --  although 
there  is  a  recent  major  thrust  in  this  direction  --  a 
source  of  lead  ultrafines  for  use  on  this  program  has 
been  difficult  to  locate. 
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Discussion  wi.th  research  leaders  in  this  field 
has  established  that  Cornell  University  has  the  capability 
of  producing  the  materials.  Cornell  may  be  the  only 
organization  in  the  United  States  capable  of  doing  this 
work,  although  similar  capability  probably  can  be  located 
in  several  foreign  research  centers.  No  negotiations 
have  been  undertaken  to  buy  the  particles  from  Cornell, 
but  preliminary  contacts  have  been  made. 

.  ;  Cryogenic  Refrigerator  Performance 

As  part  of  this  program  it  was  necessary  r,<-  v  /el  op 
a  raw  regenerator  performance  model  based  on  the  pri'  r 
work  of  several  researchers  to  predict  efficiency  vs. 
specific  heat  effects  at  low  temperature.  This  new  model, 
when  applied  to  four  test  cases  typical  of  cryogenic  re - 
f r !  go r i.  tor  performance  rear  i remap  tn ,  j>r  -'d  ’  o  * s  i  gnif  j  car  t 
improvemen  tr  in  operating  ch'.u  ar  ter  is  ti  cc  when  selected 
uitrafine  materials  are  used  to  form  the  heat  transfer 
elements  in  the  regenerator  matrix. 

In  particular,  the  model  indicates  that,  the  size  of 
the  cold  stage  mechanism  can  be  reduced,  a  significant 
safety  factor  (by  way  of  excess  refrigeration  capacity) 
can  be  provided,  a  significant  power  savings  can  be 
achieved,  or  an  optimum  combination  of  these  three  new 
degrees  of  freedom  can  be  selected.  These  savings  are 
judged  to  be  a  significant  benefit  to  existing  refrigerator 
technology . 

5 . 4  Regenerator  Element  Producibility 

Of  the  questions  addressed  during  this  research 
study,  the  one  which  is  least  confidently  answered 
involves  the  producibility  of  regenerator  elements. 

Based  on  recently  reported  measurement  results,  a  green 
compact  of  the  powdered  material  can  be  formed  without 
destroying  the  surface  enhanced  specific  heat  which  is 
characteristic  of  the  ultrafines.  The  suitability  of 
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such  a  compact  for  use  in  a  periodic  flow  regenerator, 
due  to  mechanical  strength  considerations,  can  only  be 
determined  by  experimental  simulutions  and  evaluation. 

Several  other  fabrication  techniques  have  been 
identified  and  determined  to  be  candidates  for  forming 
regenerator  elements.  Each  of  them,  however,  must  be 
further  evaluated  experimentally  because  very  little 
work  on  the  forming  and  shaping  of  sub-micron  powders 
has  been  reported  in  the  literature.  This  is  a  new 
technological  area. 

The  fabrication  techniques  selected  for  initial  use 
during  the  later  phases  of  this  program  will  be  the 
green  compact  method,  and  the  green  compact  method 
combined  with  the  mercury  porosimetry  technique  described 
in  the  text. 

5*5  Safety 

Experimental  work  with  ultrafines  of  lead  introduces 
additional  hazards  beyond  the  toxicity  problem  normally 
encountered.  The  large  surface  to  volume  ratio  of  the 
powder  results  in  a  high  energy,  exposed  surface  area. 
This  large  surface  area  will  react  readily  with  oxygen 
(or  similar  substances)  and  tends  to  spontaneously  burn 
or  explode.  This  problem  is  minimized  by  storing  and 
handling  the  materials  only  in  an  inert  gas  environment. 

The  ultrafines  pose  an  inhalation  danger,  very 
similar  to  lead  vapor,  thus  escalating  the  health  con¬ 
siderations,  and  normal  safety  precautions  require  the 
use  of  a  vented  laboratory  hood  or  glove  box,  as  well 
as  the  use  of  personal  respirators  for  workers  who 
directly  contact  the  material.  With  these  precautions, 
the  hazards  can  be  reduced  to  an  acceptably  low  level. 

5.6  Summary 

In  summary,  the  expected  small  particle  specific 
heat  enhancement  at  low  temperature  has  been  verified 
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to  exist.  New  regenerator  analysis  models  have  been 
developed  as  part  of  this  program  to  predict  the 
magnitude  of  technology  improvement  to  be  achieved 
through  use  of  these  materials  in  cryogenic  regenerators. 

It  has  been  determined  that  laboratory  techniques 
for  forming  the  particles  have  been  developed  to  an 
advanced  state  by  several  research  groups.  No  proven 
technology  currently  exists  for  forming  and  shaping 
the  particles  to  a  configuration  suitable  for  the 
cryogenic  regenerator  application,  but  no  insurmountable 
problems  in  this  development  area  have  been  identified. 

Safety  considerations  will  require  that  special 
precautions  be  observed  in  working  with  the  materials, 
but  again  no  insurmountable  problems  are  anticipated. 

Based  on  the  results  of  this  current  program  phase, 
the  achievable  improvement  in  refrigerator  performance 
which  should  be  realized  through  use  of  small  particle 
technology  warrants  a  continuation  of  this  research 
activity.  The  experiment  plan  for  later  program  phases 
is  included  as  an  appendix. 
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APPENDIX 


TEST  PLAN 

1.0  Introduction 

The  Low  Temperature  Regenerator  study  program  is 
divided  into  several  phases.  The  first  phase,  now 
completed,  encompassed  both  a  literature  survey  and  the 
development  of  analytical  models  for  small  particle 
specific  heat,  regenerator  efficiency,  and  overall  refrig¬ 
erator  performance. 

The  second  phase  is  designed  to  include  the  production, 
as  well  as  mechanical  and  thermal  testing,  of  small 
particles  and  small  particle  composites.  In  addition, 
the  regenerator  performance  efficiency  models  developed 
in  Phase  1  will  be  evaluated  and  optimized. 

The  third  phase  will  ccmpgre  computed  values  of 
improved  regenerator  performance  with  measured  results 
derived  from  a  refrigerator  test  bed  as  well  as  from 
simulated  dynamic  testing  in  a  regenerator  test  fixture 
upgraded  from  the  efficiency  model  test  cryostat  which 
will  be  used  in  Phase  2. 

This  test  plan  details  the  test  logic  and  requirements 
for  the  Phase  2  work.  Successful  completion  of  the  tasks 
outlined  herein  will  provide  the  management  and  technical 
information  necessary  to  move  into  Phase  3  and  conduct 
simulated  dynamic  testing  and  actual  refrigerator  test¬ 
ing  of  the  proposed  regenerator  design.  The  Phase  3 
program  will  evaluate  the  predictions  of  input  power 
savings  and  refrigeration  capacity  improvements  derived 
from  the  Phase  1  and  Phase  2  work. 

2.0  Scope  of  Test  Plan 

This  test  plan  addresses  verification  and  test  activities 
in  two  distinct  Phase  2  areas;  evaluation  of  small  particle 
(and  small  particle  composite)  specific  heat  enhancement, 
and  evaluation  of  the  segmented  regenerator  performance 
model  discussed  earlier  in  this  report. 
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3.0 


Technical  Requirements 


3.1  Small  Particle  Specific  Heat 

Because  of  the  lack  of  reported  measurement  results 
for  materials  which  are  candidates  for  regenerator 
matrix  packings,  the  specific  heat  of  selected  materials 
produced  in  several  batches  will  be  evaluated.  The 
sampling  strategy  will  be  designed  to  preclude  a  situation 
where  suspected  anomalies  due  to  production  peculiarities 
adversely  affect  conclusions.  Unless  toxicity,  safety, 
or  producibility  problems  force  a  reconsideration,  the 
two  test  materials  will  be  lead  and  selenium. 

3 . 2  Small  Particle  Composites 

Small  particle  composites  formed  from  two  different 
materials  and  several  different  production  techniques 
will  be  evaluated  for  mechanical  and  thermal  character¬ 
istics  at  room  temperature  and  (selectively)  at  low 
temperatures.  The  redundancy  in  this  approach  will  again 
allow  evaluation  of  any  anomalies  which  could  cloud 
results  of  the  measurements.  The  composites  will  be 
formed  as,  a)  pressed  green  compacts,  b)  mercury  impreg¬ 
nated  pressed  green  compacts,  c)  permeable  membranes 
loaded  with  the  ultrafines,  and  d)  ultrafines  dispersed 
in  a  binder  and  processed  to  wafer,  spherical  or  cylindrical 
shapes . 

3.2.1  Physical  Characteristics 

A  wafer  configuration  of  the  four  forms  of  composites 
will  be  used  for  a  preliminary  evaluation  of  physical 
characteristics  of  ultrafines.  This  evaluation  will  be 
conducted  in  two  parts.  First,  to  determine  binding 
integrity,  a  strip  of  ordinary  bond  paper  sandwiched 
between  the  weighed  wafer  and  a  100  gm  weight  will  be 
pulled  for  two  inches;  then  the  wafer  shall  be  inspected 
and  re-weighed  to  determine  loss  of  ultrafines  due  to 
abrasion.  Second,  a  fracture  test  will  be  conducted 
both  at  room  temperature  and  at  77° K.  This  test  will 
consist  of  applying  an  axial  force  (via  a  1/16"  diameter 
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ball-end  rod)  to  the  wafer,  and  incrementally  increasing 
the  load  until  significant  deformation  or  fracturing  of 
the  composite  is  observed. 


Thermal  Characteristics 

The  composites  will  be  evaluated  as  to  heat  capacity 
characteristics  in  the  range  from  5  to  100°K.  In 
addition  to  the  four  types  of  composites  defined  in 
Paragraph  3*2,  heat  capacity  measurements  will  also  be 
performed  for  up  to  two  cluster  structure  materials 
(as  defined  in  Section  II)  and/or  selected  "anomaly" 
materials  which  evidence  large  entropy  effects  in  the 
low  end  of  the  temperature  range. 


Regenerator  Efficiency  Equation 

In  order  to  test  and  verify  the  accuracy  of  the 
segmented  regenerator  efficiency  equation  developed  in 
Section  III,  a  progressive  evaluation  program  is  required. 
This  program  encompasses  the  following  four  test  sequences: 

Test  Sequence  I  -  Constant  Regenerator  Matrix 

Stream  Capacity  and  constant 
gas  stream  capacity. 

Test  Sequence  II  -  Multi-segmented  Matrix  Stream 

Capacity  and  constant  gas  stream 
capacity. 

Test  Sequence  III  -  "Continuously"  Variable  Matrix 

Stream  Capacity  and  constant  gas 
stream  capacity. 

Test  Sequence  IV  -  "Continuously"  Variable  Matrix 

Stream  Capacity  and  variable  gas 
stream  capacity. 

Sequences  I,  II,  and  the  initial  testing  of  Sequence  III 
will  be  conducted  with  the  regenerator  operating  between 
room  temperature  and  77°K.  The  remainder  of  Sequence  III 
as  well  as  Sequence  IV  will  be  conducted  with  the  regen¬ 
erator  operating  as  low  as  the  liquid  helium  temperature 
range . 
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3.3.1 


Sequence  I  Evaluation 

With  an  experimental  regenerator  (Rl)  built  up  of 

material  of  known  specific  heat»  the  Lamberts on-London 

finite  difference  method  of  calculating  efficiency,  as 

well  as  the  empirical  equation  from  Section  III  of  this 

report  will  be  evaluated.  The  reciprocity  effects  of 

and  0^  on  regenerator  efficiency  (i.e.  the  effect  on 

efficiency  of  increases  of  C  versus  decreases  of  C  ) 

r  g 

will  be  explored.  A  determination  will  be  made  whether 
or  not  axial  heat  conduction  must  be  accounted  for  in 
order  to  correlate  measured  efficiency  with  calculated 
values.  Test  conditions  will  encompass  NTU  values  up 
to  about  200  and  C ^/C  values  ranging  from  3  or  4  down 
to  below  unity. 

3*3.2  Sequence  II  Evaluation 

With  an  experimental  segmented  regenerator  (R2) 
built  up  of  two  or  three  segments  of  known  specific  heat, 
the  segmented  regenerator  model  derived  in  Section  III 
will  be  evaluated.  As  in  Sequence  I,  reciprocity  effects 
for  high  NTU  values  and  Low  ratios  will  be  invest¬ 

igated,  and  a  determination  of  the  importance  of  axial 
heat  conduction  will  be  accomplished. 

3.3.3  Sequence  III  Evaluation 

With  an  experimental  regenerator  (R3)  built  up  of  a 
larger  number  of  divided  segments,  the  segmented  regen¬ 
erator  model  of  Section  III  will  be  applied  to  approximate 
the  continuously  variable  matrix.  C  and  C  reciprocity 
characteristics  will  be  investigated  again  for  high  NTU 
value  and  low  C^Cg  ratios.  Results  will  be  compared 
with  calculated  values  from  the  segmented  regenerator 
model,  and  the  importance  of  axial  heat  conduction  on 
performance  predictions  will  be  assessed. 

3.3.4  Sequence  IV  Evaluation 

The  regenerators  used  in  the  Sequence  III  program 
will  be  re-evaluated  in  the  77  -  5° K  range  where  the 
specific  heat  and  thus  the  stream  capacity  of  helium 
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gas  is  a  function  of  temperature  at  the  pressure  of 
interest.  As  before,  test  results  will  be  compared 
with  calculated  predictions,  and  the  importance  of  axial 
heat  conduction  will  be  evaluated.  During  a  sub-part  of 
this  sequence  the  helium  gas  flow  rate  will  be  varied  in 
a  linear  manner  so  that  a  preliminary  model  for  predict¬ 
ing  time  varying  flow  and  phase  angle  effects  can  be 
empirically  compiled. 


Test  Sequence  Requirements 


Specific  Heat  Measurements 

The  planning  network  for  specific  heat  measurements 
is  shown  in  Figure  32.  The  schedule  of  tests  depicted 
therein  is  designed  to  provide  the  information  base 
necessary  to,  a)  verify  the  magnitude  of  surface  enhanced 
specific  heat  for  powders  and  powder  composites,  and  b) 
perform  a  preliminary  evaluation  of  the  mechanical 
characteristics  of  powder  composites  for  use  in  regener¬ 
ator  fabrication.  This  information  will  be  used  to 
define  the  scope  of  follow-on  Phase  3  work  and  to 
evaluate  whether  to  design  and  build  a  prototype 
regenerator  using  these  effects. 


Description  of  Test  Activity  Blocks  (per  Figure  32) 


Block  No. 

Title 

Discussion 

1 

Procure  Powders 

Ultrafine  powders  meeting 
the  size  distribution  re¬ 
quirements  of  Section  II  of 
this  report  will  be  procured. 

2 


Cellulose  acetate -based 
and  teflon-based  membrane 
materials  with  defined  pore 
sizes  will  be  procured. 


Procure  Basic 
Membrane  s 


Block  No.  Title  Discussion 

3  Build  Specific  A  specific  heat  cryostat  will 

Heat  Cryostat  be  fabricated  from  the  design 

package  developed  as  part  of 
the  Phase  I  program. 

4  Calibrate  Cryostat  Using  a  specimen  of  known 

specific  heat,  the  cryostat 
will  be  calibrated  in  the 
5  -  100 °K  range. 

5  Preliminary  Mem-  The  pressure  drop  of  gaseous 

4  » 

brane  Evaluation  He  a)  through  a  10  cm 

cylindrical  stack  of  membrane 
wafers  and  b)  through  a  wound 
spiral  cylinder  10  cm  in  lengh 
will  be  determined.  Coarse 
determination  of  brittleness/ 
durability  at  78°K  will  be 
performed. 

6  Decision  Block  Based  on  pressure  drop  and 

durability  tests  a  go-no  go 
decision  will  be  made  with 
regard  to  further  fabrication 
and  testing  of  membrane 
composites. 

7  Size  Distribution  After  receipt  of  ultrafines 

Assessment  from  the  vendor  and  just 

prior  to  further  processing 
or  use,  a  size  distribution 
measurement  shall  be  performed. 
Results  will  be  compared  with 
vendor  data  to  assess  long 
term  agglomeration  character¬ 
istics. 
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Discussion 

Using  the  calibrated  cryostat 
the  specific  heat  of  several 
batches  of  the  various  ultra- 
fine  materials  shall  be 
determined.  Statistical 
design  shall  be  employed  to 
determine  the  number  of  sets 
of  measurements  required. 

Press  Compacts  Several  compacts,  the  actual 

quantities  determined  by  stati¬ 
stical  analysis,  shall  be 
pressed  from  batches  of  each 
material  being  evaluated. 

The  binding  integrity  and 
fracture  tests  of  Para.  3*2.1 
will  be  conducted  both  before 
and  after  cold  cycling  the 
composites  to  liquid  helium 
temperature . 

Decision  Block  Based  on  an  assessment  of  the 

composite's  mechanical  charact- 
erstics,  a  decision  will  be 
made  whether  or  not  to  continue 
to  evaluate  the  material  for 
use  in  regenerator  matrix 
fabrication. 

Several  composites,  the  actual 
quantity  determined  by  statist¬ 
ical  analysis,  shall  be  formed 
in  wafer  or  cylindrical  shapes 
by  dispersing  ultrafine  powders 
in  a  binder  solution. 


Form  "Binder" 
Composites 


Mechanical 

Characteristics 


Title 

Specific  Heat 

Measurement 

(Powder) 


Block  No.  Title 


Discussion 

13  Mechanical  The  binding  integrity, 

Characteristics  particle  adhesion  character¬ 

istics,  and  formability  will  be 
evaluated  both  before  and  after 
cold-cycling  the  composite. 

14  Decision  Block  Based  on  an  assessment  of  the 

composite's  mechanical  character¬ 
istics,  a  decision  will  be  made 
whether  or  not  to  continue  to 
evaluate  the  material  for  use 
in  regenerator  matrix  fabri¬ 
cation. 

15  Mercury  Using  mercury  porosimetry 

Impregnation  techniques,  a  selected  number  of 

pressed  compacts  will  be  impreg¬ 
nated  with  mercury  so  as  to 
completely  fill  the  void  volume. 

16  Mechanical  The  binding  integrity  and 

Characteristics  fracture  tests  of  Para.  3.2.1 

will  be  conducted  both  before 
and  after  cold  cycling  the 
composites  to  liquid  helium 
temperature . 

17  Decision  Block  Based  on  an  assessment  of  the 

composite's  mechanical  character¬ 
istics,  a  decision  will  be  made 
whether  or  not  to  continue  to 
evaluate  the  material  for  use  in 
regenerator  matrix  fabrication. 

18  Produce  Loaded  Membrane  materials  will  be  loaded 

Membranes  with  ultrafine  powders.  For 

cellulose-based  membranes,  the 


Block  No. 
18 

continued 


19 


20 


21 


22 


Title 

Produce  Loaded 
Membranes 


Mechanical 


Decision  Block 


Specific  Heat 
Measurement 
(Cluster 
Materials) 


Specific  Heat 
Measurement 
( "Anomaly" 
Material) 


Discussion 

forming  solution  will  be  loaded 
to  a  high  density  prior  to  curing 
as  a  membrane.  The  teflon-based 
membrane  material  will  be  loaded 
by  fusing  techniques  after 
formation,  unless  a  more  efficient 
process  can  be  developed  by  the 
supplier. 

The  binding  integrity,  particle 
adhesion  characteristics,  and 
formability  will  be  evaluated  both 
before  and  after  cold-cycling  the 
composite . 

Based  on  an  assessment  of  the 
composite’s  mechanical  character¬ 
istics,  a  decision  will  be  made 
whether  or  not  to  continue  to 
evaluate  the  material  for  use  in 
regenerator  matrix  fabrication. 

Using  the  calibrated  cryostat, 
the  specific  heat  of  selected 
cluster  materials  (see  Section  II) 
in  both  bulk  and  wire  screen  form, 
if  available,  will  be  measured. 

Using  the  calibrated  cryostat, 
the  specific  heat  of  selected 
materials  showing  an  entropy 
anomaly  in  the  temperature  range 
of  interest  will  be  determined. 
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Block  No.  Title  Discussion 

23  Size  Distribution  A  size  distribution  determination 

Assessment  of  samples  of  ultrafines  from 

several  batches  of  each  material 
shall  be  performed.  Results 
shall  be  compared  with  the  earlier 
measurements  (Block  7)  to  further 
assess  agglomeration  characteristics. 


24 


25 


26 


Analysis  of  An  analysis  shall  be  performed 

Powder  Specific  to  determine  if  the  achieved 
Heat  surface  enhanced  specific  heat 

correlates  with  predicted  values 
as  determined  from  the  Section  II 
equation. 


Specific  Heat 
Measurement 
(Composite 
Materials) 


Using  the  calibrated  cryostat, 
the  specific  heat  of  samples  of 
each  composite  material  shall  be 
determined.  Statistical  design 
shall  be  employed  to  determine 
the  number  of  sets  of  measure¬ 
ments  required. 


Selection  of  Based  on  the  measurements  and 

Regenerator  observations  collected,  the  best 

Materials  candidate  materials  for  use  in 

regenerator  fabrication  will  be 
selected. 


4.2  Regenerator  Efficiency  Measurements 

The  planning  network  for  regenerator  efficiency  model 
evaluation  is  shown  as  Figure  3 3.  The  schedule  of  tests 
depicted  therein  is  designed  to  provide  the  information 
base  necessary  to,  a)  verify  the  validity  of  the  segmented 
regenerator  model  (see  Section  III),  b)  evaluate  the 
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importance  of  axial  heat  conduction  and,  c)  investigate 
the  reciprocity  characteristics  of  Cr  and  Cg  variations. 
This  information  will  be  used  to  verify  the  predictions 
of  improved  refrigerator  performance  derived  from  the 
small  particle  enhanced  specific  heat. 


4.2.1 


Description  of  Test  Activity  Blocks  (per  Figure 


Block  No.  Title 


Discussion 


Measure  Test  Using  the  calibrated  specific 

Material  Specific  heat  cryostat,  the  heat  capacity 
Heat  of  screen  and/or  sphere  materials 


to  be  used  in  this  test  program 
will  be  evaluated  in  the  5°K  to 
77°K  range. 


2  Build  Regenerator  Based  on  the  design  package 

Test  Cryostat  developed  during  Phase  1,  the 

regenerator  test  fixture  will 
be  constructed. 


3  Fabricate  A  regenerator  body  will  be  fabri- 

Regenerator  (Rl)  cated  and  150,  200  and  325  mesh 

screens  of  two  different  test 
materials  (stainless  steel  and 
phosphor  bronze)  will  be  punched 
to  size ,  so  that  the  matrix 
packings  can  be  conveniently 
interchanged. 


4  High  Temperature  The  regenerator  test  cryostat 

Baseline  Test  will  be  baselined  in  the  300  -  77°K 

range  using  the  Rl  regenerator 
housing  with  .005  -  .009”  dia. 
lead  spheres,  as  well  as  150 
mesh  screen  packing.  Test 
results  will  be  compared  where 
possible  with  similar  measure¬ 
ments  by  Gifford^  and  Lechner^2^ 
in  order  to  assess  the  fixture 
performance. 
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Block  No.  Title 


Measure  Regener¬ 
ator  R1 

Efficiency  (Test 
Sequence  I) 


Low  Temperature 
Baseline  Test 


Fabricate 
Regenerator  R2 


Fabricate 
Regenerator  R3 


Discussion 

The  R1  regenerator  using  all 
combinations  of  screen  mesh 
materials  and  sizes  will  be 

evaluated  in  the  300  -  77°K 
range ,  Square  wave  flow  patterns 
will  be  utilized  and  peak  helium 
flow  rates  will  be  varied  in  the 
2.5  -  25  SCFTC  range.  Blow  times 
will  be  varied  from  50  -  500  ms. 
Efficiency  must  be  determined 
to  three  significant  figures. 

The  regenerator  test  cryostat 
will  be  baselined  in  the  77  -  5°K 
range  using  the  R1  regenerator 
housing  with  .005  -  .009"  dia. 
lead  sphere  packing.  Test  results 
will  be  evaluated  against  the 
calculated  values  extrapolated 
from  Block  4  results, 

A  regenerator  will  be  fabricated 
in  such  a  manner  that  it  can  be 
partitioned  into  two  (as  well 
as  three)  segments  of  matrix 
material.  Each  segment  will  be 
packed  with  a  matrix  of  pre¬ 
selected  heat  oapacity. 

A  regenerator  will  be  fabricated 
in  such  a  manner  that  it  can  be 
partitioned  into  more  than  three 
sepients  of  matrix  material. 

Eaoh  segment,  not  necessarily 
of  equal  lengths,  will  be  packed 
with  a  matrix  of  pre-selected 
total  heat  capacity, 
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Title 


Discussion 


Verify 
Lambertson- 
London  Model 


Measure  Regen¬ 
erator  R2 
Efficiency  (Test 
Sequence  II) 


Measure  Regen¬ 
erator  R3 
Efficiency  (Test 
Sequence  III) 


Verify  Segmented 

Regenerator 

Model 


Using  the  data  obtained  from 
Block  5»  the  C/C  capacity 

^  O 

ratio  effect  on  efficiency  will 
be  investigated.  Effects  of 
axial  heat  conduction  and  wall 
effects  will  be  considered. 

The  R2  regenerator  using  appro¬ 
priate  combinations  of  screen 
mesh  materials  and  sizes  will  be 
evaluated  in  the  300  -  77° K 
range.  Square  wave  flow  patterns 
will  be  varied  in  the  2.5  -  25 
SCFM  range.  Blow  times  will  be 
varied  from  50  -  500  ms. 

Efficiency  must  be  determined 
to  three  significant  figures. 

The  R3  regenerator  using  appro¬ 
priate  combinations  of  screen 
mesh  materials  will  be  evaluated 
in  the  300  -  77°K  range.  Square 
wave  flow  patterns  will  be 
utilized  and  peak  helium  flow 
rate  will  be  varied  in  the 
2.5  -  25  SCFM  range.  Blow  times 
will  be  varied  from  50  -  500  ms. 
Efficiency  must  be  determined  to 
three  significant  figures. 

Using  the  data  obtained  from 
Blocks  10  and  11,  the  C^Cg 
capacity  ratio  effect  on  segment 
efficiency  and  overall  regenerator 
performance  will  be  investigated. 
Effects  of  axial  heat  conduction 
will  be  considered. 


Block  No. 

Title 

13 

Measure  Low 

Temperature 

Regenerator 

Efficiency  (Test 

Sequence  IV ) 

Discussion 

Using  preselected  combinations 
of  matrix  packings  in  the  R2  and 
R3  regenerators,  performance  will 
be  evaluated  in  the  77  -  5°K 
range.  Square  wave  flow  patterns 
will  be  utilized  and  peak  helium 
flow  rates  will  be  varied  in  the 
2.5  -  25  SCFM  range.  Blow  times 
will  be  varied  from  50  -  500  ms. 
Efficiency  must  be  determined  to 
three  significant  figures. 


14  Final  Model  Using  the  data  obtained  from 

Assessment  Block  13,  the  Cj/Cg  capacity 

ratio  effect  on  low  temperature 
regenerator  efficiency  and  overall 
performance  will  be  investigated. 
Effects  of  axial  heat  conduction 
will  be  considered. 


(1)  W.  E.  Gifford,  Basic  Investigation  of  Cryogenic 
Refrigeration  Methods,  Air  Force  Flight  Dynamics 
Laboratory,  Wright  Patterson  Air  Force  Base,  Ohio, 
45433,  AFFDL-TR-68-61  (June  1968). 

(2)  R.  A.  Lechner,  Investigation  of  Regenerators  and 
Pulse  Tube  Cryogenic  Coolers,  United  States  Army 
Electronics  Command,  Fort  Monmouth,  N.J., 

EC0M-3409  (May  1971). 
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